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ABSTRACT 


In satellite laser ranging, atmospheric refraction increases the optical 
pathlength from the ground to an orbiting satellite. Two-color laser ranging 
systems can be used to determine the atmospheric delay by measuring the 
difference in propagation times between optical pulses transmitted at the two 
different wavelengths. The differential propagation time can also be used to 
infer surface pressure. 

The statistical properties of the signals reflected from the 
re :ro-ref lector equipped satellites are studied. It is found that coherence 
interference between pulse reflections from retro-reflectors of different 
ranges on the array platform i6 the primary cause of signal fluctuations. 

The performance of a cross-correlation technique for estimating the 
differential propagation time is analyzed by considering both shot noise and 
speckle. For the retro-reflector arrays, timing performance is dominated by 
interference induced speckle, and the differential propagation time cannot be 
resolved to better than the pulse widths of the received signals. , 

The differential timing measurements obtained over a horizontal path are 
analyzed. The results are used to verify both the two-color ranging 
measurement technique and the cross-correlation timing algorithm. 

The ocean-reflected pulse measurements obtained from the airborne 
two-color laser altimeter experiment are presented. The reflected pulse shapes 
are used to infer the 6ea state and the ocean wind speed. The measured 
differential propagation times are used to verify the pressure measurement 


technique . 
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1 . INTRODUCTION 

The feasibility of using pulsed lasers to range to earth-orbiting 
satellites was first demonstrated by researchers at the NASA-Goddard Space 
Flight Center in 1964 (Plotkin et al., 1965 J. Since then, this technique 
has been used by numerous groups for applications in global surveying (Smith 
et al., 1972J, geophysical research [Smith et al , 1973], and positioning 
(Smith et al., 1973]. In these applications, pulsed laser ranging systems 
are used to estimate the distances from ground-based ranging sites to retro- 
reflector equipped satellites by measuring the roundtrip optical propagation 
times. The accuracies of these systems are limited by the effects of 
atmospheric refraction, which increases the average optical pathlength be- 
tween the ranging site and the orbiting satellite by 2.5 meters when the 
satellite is near zenith and by more than 13 meters when it is at 10° eleva- 
tion [Abshire and Gardner, 1985]. Several formulas have been developed 
which can correct range errors to centimeter level for the single-color 
systems [Saastamoinen, 1972; Marini and Murray, 1973; Gardner, 1977; Davis 
et al., 1984; Treuhaft et al., 1984]. However, these formulas require 
accurate meteorological measurements at the ranging site during the satellite 
pass, which may not be feasible, particularly in certain inaccessible regions of 
the earth. An alternate approach is to employ two-color laser ranging systems 
to determine the atmospheric correction directly by measuring the propagation 
time difference induced by dispersion between pulses at two different optical 
frequencies (Prilepin, 1957; Bender and Owens, 1965]. Because the atmospheric 
dispersion in the optical frequency band is small, the differen xal propagation 
times must be measured with picosecond-level accuracies if cm-level ranging 
accuracies are desired [Bender and Owens, 1965]. 
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More recently, the multicolor laser ranging technique has been extended 
to applications in remote sensing of atmospheric pressure. By using the 
facts that the atmosphere is dispersive, and the group refractivity is pro- 
portional to the surface pressure, Gardner [1979] proposed a remote sensing 
technique to globally measure atmospheric pressure using a two-color laser 
altimeter. Making pressure measurements over the ocean using this remote 
technique is highly desirable because the weather over land usually origina- 
tes from oceans. Accurate measurements of ocean surface pressure would 
allow better seasonal climate modelling and long-range weather forecasting. 

In general, when the altimeter uses 1.064 - and 0.353-pm YAG laser wave- 
lengths, picosecond timing accuracies are required for millibar-level 
pressure accuracies. 

The differential timing accuracy, in general, is strongly dependent on 
the orientation and the local height distribution of the target surface pro- 
file with respect to the optical axis. Experimental results obtained using 
high speed detectors indicate that picosecond timing accuracy is achievable 
for specular targets like non-dispersive retroref lectors [Abshire and 
Kalshavon, 1983; Im et al., 1983; Abshire and Gardner, 1985], However, these 
target geometries constitute only a small fraction of the targets which are 
available. In order to extend the applicability of the two-color ranging 
technique, targets with randomly varying local surface profiles, such as the 
ocean and satellites containing retroref lector arrays, must also be con- 
sidered. The reflections from such targets present three major difficulties 
for timing estimation. First, the reflected pulses are broadened due to the 
range spreads of the target surface profiles. Second, the time-resolved 
speckle fluctuations within the received pulses will distort the pulse shapes 
when the range spreads of the targets are longer than the transmitted 
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laser pulse widths (Tsai and Gardner, 1985]. Third, the mean received pulse 
shapes from randomly changing target surfaces are generally not known 
a priori so that conventional optimal tim: ng estimators cannot be used. 

The achievement of picosecond timing accuracy when ranging to targets of 
opportunity is, therefore, a challenging task which requires careful 
investigation of both the reflection process and the ranging receiver 
design. This dissertation is centered around these two aspects, with empha- 
sis on the two-color ranging accuracy in the presence of major noise sources 
such as shot noise, speckle, and background noise. 

Chapter 2 reviews the established theories behind the pulsed laser 
ranging technique, the two-color atmospheric correction technique, and the 
two-color pressure measurement technique, and their corresponding accura- 
cies. 

In the first portion of Chapter 3 the statistics of the diffuse target 
reflected laser signal (Gardner, 1982] is summarized. The statistics of the 
signals reflected from the cube-corner reflector arrays is studied, using 
partially developed speckle theory and Fresnel diffraction. This study has 
practical application in ranging to the existing retroref lector-equipped 
satellites such as the Laser Geodynamic Satellite (LAGEOS) and Geodetic 
Earth Orbiting Satellites (GEOS). 

In addition to the range errors induced by the atmosphere, the perfor- 
mance of the laser ranging systems is also limited by the presence of noise. 
In Chapter 4 the performance of the Maximum-Likelihood timing estimators are 
investigated under the influence of the time-resolved speckle and background 
radiation. 

When the targets of interest have temporally varying local surface 
profiles, it is not possible to predict the mean reflected signal a priori. 
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Consequently, the Maximum Likelihood estimators, which require knowledge of 
the mean signal, cannot be used. In Chapter 5 a two-color cross-correlation 
estimator is studied and its corresponding performance is analyzed. The 
results are applied to the pulse reflections from ground diffuse targets, 
earth-orbiting CCR arrays, and the ocean. 

For high precision two-color laser ranging, the transmitted laser beams 
at the two frequencies must be aligned so that the reflected signals are 
highly correlated. Any misalignment will cause decorrelation of the laser 
pulse shape and degrade the timing performance of the cross-correlation esti- 
mator. In addition, the relative displacement between the transmitter and 
receiver axes can influence the amount the signal strength received by the 
system, which in turn affects the timing performance because the timing 
accuracy is directly proportional to the signal strength. In Chapter 6 we 
investigate the system alignment effects on the performance of the cross- 
correlation differential timing estimator. 

In parallel with the theoretical research on the two-color cross- 
correlation technique, experimental research has been undertaken at NASA- 
Goddard Space Flight Center to verify the feasibility of this technique under 
actual field conditions. In Chapter 7 we report on the recent two-color range 
measurements made over the horizontal paths. These measured results are ana- 
lyzed and compared with the predictions based upon the theory of Chapter 5. 

An airborne two-color altimeter experiment was conducted in April of 
1985 at NASA/Goddard-Wal lops Flight Facility in which the ocean-reflected signal 
measurements were collected. The primary objectives of this mission were to 
verify the two-color pressure measurement technique and the cross- 
correlation technique. The processed results are analyzed and discussed in 


Chapter 8. 
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2. RANGING AND PRESSURE MEASUREMENT TECHNIQUES : REVIEW 

2 . 1 Introduction 

Both single- and two-color satellite laser ranging (SLR) systems can be 
used to estimate distance by measuring the round-trip propagation times of the 
laser pulses. However, the performance of such systems is limited by the 
effects of atmospheric refraction, which increases the optical pathlengths 
between the satellite and the ranging site, and the effects of atmospheric 
turbulence, which introduces random fluctuations in the refractivity along the 
propagation path. Atmospheric refraction increases the average round-trip 
optical pathlength to an orbiting satellite by 2.5 m when it is near zenith and 
by more than 13 m when the satellite is at 10* elevation. Strong turbulence 
can introduce random errors of about several cm at 10’ elevation [Abshire and 
Gardner, 1985]. Consequently, atmospheric refraction is the most dominant 
source of error in satellite ranging. Normally, the error contributions due 
to turbulence are negligible. 

2.2 Ranging Geometry 

The geometry for a satellite laser ranging system and its ground target 
is shown in Fig. 2.1. The round-trip straight-line pathlength R g is the 
parameter of interest. However, due to the atmospheric refraction effect, the 
round-trip optical pathlength R is longer than R . The atmospheric correc- 
tion AC (that is, the range error due to atmospheric refraction) is the dif- 
ference between the round-trip optical pathlength and the straight-line 
pathlength : 

AC - (R - R ) . 

O 6 


( 2 . 1 ) 


Satellite 
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The optical pathlength between the satellite and the ground target is defined 
as the integral of the group refractive index along the ray path: 

R - 2 / n dl (2.2) 

o J c g 

where dl is an incremental length along the ray path c, and n^ is ths group 
refractive index. This refractive index is related to the group refractivity 
by [Marini and Murray, 1973] 

n - 1 ♦ l(f 6 N . (2.3) 

g g 

If the horizontal refractivity gradients are small, the ray path will lie 
entirely in a plane. In this case, the optical pathlength for a pulsed laser 
system is given by [Marini and Murray, 1973] 


R 

o 


2 J 


r (1 ♦ 10" 6 N ) 

sat g 


8 in 6 


dr 


earth 


(2.4) 


where r is the geocentric altitude, and 6 is given by Snell's Law for a 
spherical stratified medium. The atmospheric correction in the spherically 
symmetrical atmosphere is thus 


AC 


2 I 


sat 

earth 



s in 


N 

e 


dr 


♦ 



sat 

earth 




(2.3) 


The first term in Eq. (2.5) is the velocity correction error, while the second 
term is the difference bctweer the geometric length of the ray and the 
straight-line path. 
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2 . 3 Atm o spheric Correction Model For Single-color Laser Ranging 

The atmospheric correction can be evaluated by using an appropriate model 
for the group ref ractivity. For typical refractivity profile that includes 
the nonsymmetric behavior of tnc atmosphere, AC can be modeled approximately 
as [Gardner and Axford, 1980 J 

8 , 8 ? 8 -, cos a 8 , sin a 

AC = — + — =y~ + - + (2.6) 

sin E sin J E sin E tan E sin E tan E 

where E is the elevation angle, a is the azimuth angle, and the coefficients 
3 ^, 8 -,. 83 * 84 are functions of the refractivity profile and laser wavelength. 
The first two terms in Eq. (2.6) are spherical correction terms corresponding 
to a spherically symmetric atmosphere, while the last two terms include the 
effects of horizontal refractivity gradients. Marini and Murray 1 1 973 J have 
developed accurate formulas for calculating 8 ^ and in terms of the surface 
pressure, temperature and water vapor pressure at the ranging site during the 
satellite pass. More recently, Gardner [1977] has derived expressions for 83 
and 8 ^ in terms of the horizontal pressure and temperature gradients at the 
ranging site. Both formulas have been checked extensively by comparing them 
with the ray tracing data. The results indicate that the error in approxi- 
mating the atmospheric correction with a formula of the form given by 
Eq. (2.6) is less than a few centimeters for elevation angles above 20°. 

If ranging accuracies of a few centimeters or less are desired, an 
extensive weather station network may be required to obtain sufficiently 
accurate measurements of the meteorological parameters during the satellite 
pass [Gardner et al., 1978J. This requirement may limit the applicability of 
che single-color laser ranging systems. 
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2 .4 Two-color Laser Ranging Technique 

Two-color laser ranging provides an attractive alternative to models for 
determining the atmocpheric correction because requirements on surface 
meteorological data can be eliminated. The technique is based upon the fact 

* 

that the dispersion in the atmosphere is frequency dependent. This causes the 
optical pathlength at the two different frequencies to be proportional to the 
path-integrated atmospheric density. Therefore, the difference in optical 
pathlengths at two laser frequencies is a measure of the refractive conditions 
existing over the propagation path at the instant the measurements are taken 
and can be used to estimate the atmospheric correction [Bender and Owens, 

1965]. 

Let R q and R q denote the round-trip optical pathlengths at wavelengths X^ 
and X^ respectively. Then the atmospheric correction at wavelength Xj is 
approximately 

AC. = Y (R - R ) = Y AR (2.7) 

2 1 

where 


Y = (n„ - 1) / (n - n ) 


g 


1 


g* 


g 


( 2 . 8 ) 


and n and n are £he group refractive indices of the atmosphere measured at 
8 1 g 2 

wavelengths X^ and , respectively. Equation (2.7) is accurate at optical 
frequencies. If we ignore the small water vapor effects at optical 
frequencies, Y can be written as 


Y = f(X 2 ) / (f(X 2 ) - f ( Xj )) 


(2.9) 


where 
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f(A.) = 0.965 + 0.0164 2 ♦ 0.000228 A." 4 , i*l,2, (2.10) 

and A^ is expressed in units of micrometers. The function f(A^) accounts for 
the atmospheric disperion which is responsible for the difference in oath- 
lengths at the two frequencies. It is derived by curve fitting the actual 
experimental measurements of the dispersion. In Table 2.1 y is listed for all 
pairwise combinations of the fundamental (1.064-pm), doubled (0.532-lim) and 
tripled (0.355-ym) YAG laser wavelengths. 

The RMS error in atmospheric correction (o _ ) is related to the RMS 
error in differential pathlength (o^) through 


AC, 


■ Y o 


AR 


( 2 . 11 ) 


Because Y is on the order of ten for all wavelength pairs, it is necessary to 
determine the differential pathlength with an accuracy which is approximately 
ten times greater than the desired accuracy for the atmospheric correction. 
Therefore, in order to achieve cm-level accuracy in atmospheric correction, 
which is typical by using a single-color system, the differential pathlength 
must be measured to an accuracy of a few millimeters, or equivalently, the 
differential propagation time of the optical pulses at the two frequencies 
must be measured to an accuracy of a few picoseconds. 

This accuracy requirement can be eased considerably if the results of 
many measurements during the satellite pass are averaged. However, the 
satellite position must be taken into account when measurements are averaged 
because both AC^ and AR are functions of the satellite azimuth and elevation 
angles. In general, regression models are used to fit the calculated values 
of the atmospheric correction to Eq. (2.6) before averaging [Gardner and 
Axford , 1980]. 
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TABLE 2.1. 

ATMOSPHERIC CORRECTION TO DIFFERENTIAL PATHLENGTH RATIO 



X 2 (yn) 

Y 

1.064 

0.532 

21.1 

1.064 

0.355 

7.45 

0.532 

1.064 

-22.1 

0.532 

0.355 

12.1 

0.355 

1.064 

- 8.45 

0.355 

0.532 

-13.1 
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2.5 Surface Pressure Measurement Technique 

Although the atmospheric correction depends on the atmospheric pressure, 
temperature, and humidity along the propagation path, theoretical results show 
that it is most sensitive to pressure [Gardner, 1979]. T n fact, when the 
laser is pointed at nadir, AC is approximately proportional to the atmospheric 
pressure at the earth's surface. Therefore, surface pressure can be computed 
from measurements of AC. 

If the straight-line pathlength (R ) between the satellite and the laser 

s 

footprint on the earth's surface is known, AC can be calculated by using the 
single-color timing measurements and Eq. (2.4). Because R must be known to 

O 

within a few centimeters for the approach to be effective, it is probably not 
practical. An alternate approach is to use a two-color laser ranging system 
to determine AC by calculating the differential pathlength measured at the two 
frequencies : 

AR » AC X - AC 2 . (2.12) 

In general, Eq. (2.6) must be used to determine AK. However, for 
pressure measurement applications we are only interested in situations where 
the nadir angles are small. In this case, the gradient correction terms given 
by Eq. (2.6) contribute less than 1 cm as compared with 2- to 5-meter 
contribution from spherical connection terms. Consequently, the more 
simplified model given by Eq. (2.5) can be used for pressure measurement 
applications without loss of accuracy. If Eqs. (2.5) and (2.12) are evaluated 
and then solved for the atmospheric pressure at the ranging site, we obtain 


[Gardner, 1979] 
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2 L/o 

v m -b + (b - 4cd) L 

2a 


(2.13) 


a 


4.73 x 10 


f(X 1 ) - f(x 2 ) 2 

T sin 2 E 3 - K 


(2.14) 


F(0 c> H) AR sin E 
2[f(X 1 ) - f(X 2 )J 


2.24 x 10 _A e 


(2.15) 


o K 7 ,n"3 , 1.084 x 10" 8 TK 
-2.357 x 10 + r 

tan E 


-13 2 2 

1.5 x 10 0 T K 

. 2-K 

sin E 


(2.16) 


F(6 c , H) * 1 + 0.0026 cos (29 c ) - 0.0003 (H) . (2.17) 

K =» 1.163 + 0.00968 cos (26^ - 0.00104 T + 0.00001435 P , (2.18) 

where 

e = water vapor pressure at the laser footprint (mbar), 

P = atmospheric pressure at the laser footprint (mbar), 

T = temperature at the laser footprint (°K), 

6 c • co-latitude cf the laser footprint, 

H * altitude of the laser footprint above sea level (km), 
and the remaining parameters were described in earlier sections of this 
chapter. 

2.6 Pressure Measurement Error 

In a realistic system it will not be possible to determine the parameters 
AR, E, e, and T exactly. Therefore, it is important to assess the effect of 
errors in these parameters on the accuracy of the recovered pressure. An 
approximate expression for the RMS pressure error is [Gardner, 1979] 
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(2.19) 


where o^, 0£. o e , 0^ are the RMS errors in AR, E, e, and T, respectively. 
The partial derivatives can be calculated by taking the derivatives of Eq . 
(2.13). The results are 


_3P 

3AR 


0.212 sin E 
f(X 1 ) - f(A 2 ) 


(mbar/mm) 


( 2 . 20 ) 


3P 

3E 


— — ^ (rabar/mrad) 
tan E 


( 2 . 21 ) 


3 0.093 (mbar/mbar) 


( 2 . 22 ) 


3P 

— 2 
3T 


2.55 x 10 _6 P 

► 2 _ 
tan E 


3.4 x 10 _5 P 2 4.2 x 10 -11 PT . . . 

+ + (mbar/ c) 


2 . 2 
T sin^E 


sin^E 


(2.23) 


Although the pressure sensitivity to temperature error is greatest at the 
lower elevation angles, its value above 20* is quite small. At 20° elevation, 
a temperature error of 50°C would contribute less than 1 mbar to the pressure 
error. Thus, only a crude estimate of the surface temperature is required. 

An accuracy of 20-30*C is easy to obtain and should be adequate. 

The pressure sensitivity to errors in water vapor pressure is constant 
with respect to elevation angle. A 10-mbar error in water vapor pressure will 
contribute approximately 1 mbar to the surface pressure error. Since water 
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vapor pressure can approach 40 rabar when the surface temperature and relative 
humidity are high, its effects cannot be ignored. Fortunately, with the 
development of the microwave radiometry, sufficiently accurate information on 
water vapor pressure is accessible. Current radiometers can measure this 
parameter to a 90% accuracy [Prabhakara et al., 1982], For an extremely 
humid condition of 40 mbar, a 10% error in water vapor pressure measurement 
corresponds to about 0.4 mbar error in surface pressure. Consequently, by 
using information collected by radiometers, the error due to water vapor 
pressure can be reduced to an insignificant amount. 

The pressure sensitivity in elevation angle is only significant at the 
lower elevation angles. At 40* or above, however, a 1-mrad error in elevation 
angle contributes less than 1 mbar error in surface pressure. Since pointing 
accuracies of 100 prad or less are easily obtained with current technology 
[ Fitzmaurice , 1975], the resulting pressure errors should be less than a few 
tenths of a millibar. 

The dominant error source of this technique is the differential path- 
length measurement. The sensitivity depends both on the choice of wavelengths 
and on the elevation angle. At low elevation angles the laser pulses travel 
through more of the atmosphere so the pressure effects are more significant. 

In addition, AR is greater for larger wavelength differences because of 
greater differences in group refractivity . 3P/8AR is plotted vs. elevation 
angle in Fig. 2.2 for three possible wavelength combinations of fundamental, 
doubled, and tripled ND:YAG laser frequencies. When ranging at nadir using 
the 0.355- and 1.064-pm wavelengths, the required differential pathlengch 
accuracy is 0.6 mm/mbar, which requires a timing accuracy of 2 psec/mbar. 


\ 
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Fig. 2.2. 


Pressure measurement sensitivity to differential pathlength 
errors as a function of the satellite elevation angle. 
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Although the pressure sensitivity to errors in AR decreases as the 
elevation angle is decreased, other error factors become more significant. 
Consequently, the best performance will be obtained with the laser pointed at 
nadir. When the laser is pointed to within a few degrees of nadir, Eq. (2.19) 
can be approximated by 


o p (mbar) 


0.212 sin E 



(mm) 


(2.24) 



• • 
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3. STATISTICS OF THE DETECTED SIGNAL 


3 . 1 Introduction 

The accuracy of the differential propagation time measurement depends 
upon the temporal characteristics of the reflected pulses and the timing 
algorithms employed in the receiver. The temporal characteristics of the 
reflected pulses are closely related to the statistics of the speckle 
fields, which are generated after the transmitted signals are reflected by the 
targets of interest. Due to the nature of the target surface profiles, the 
resultant speckle can either be fully developed or partially developed. 
Consequently, the corresponding detected signals will have different sta- 
tistics with respect to different target surface profiles. In this chapter, 
both the fully and partially developed speckle statistics will be discussed. 

3.2 Statistics of the Reflected Signals from the Diffuse Targets 

When the target surface is rough on the optical scale, and when 

sufficiently large numbers of independent scattering areas exist within 
the laser footprint, then the reflected speckle field is usually regarded as 
fully developed and follows circular complex Gaussian statistics [Goodman, 

1965 and 1975J. Gardner [ 1 982 J and Tsai and Gardner [1982] have investigated 
the temporal characteristics of optical pulses that were reflected from both 
the diffuse ground targets and the ocean. For normal incidence of the laser 
beam on the target and direct detection of the reflected pulses, the con- 
ditional mean detected signal (<S(t)>) and the auLocovariance function 
(Cg(t^,r^)> are given respectively by 

<S(t)> - <N> / d 2 p b 2 (p,z) |f(t - *(P))| 2 * h(t) , (3.1) 


and 
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C S ^ C 1 » c 2 ) * <N> ^ 2 — b ?^» z ^ / dx|f(T - <Kjj))| 2 h( t j - t) h( t. 


♦ <N> 2 K 1 / d 2 _p b 4 (_p,z) g(tj - <{/(j>)) g(t 2 - <K_p)) 


T) 

(3.2) 


where 


b n (_p,z) - |a(_p,z)| n B^^tj)) / / d 2 £ |a(_p,z) | n 


(3.3) 


♦(.£) 


2z + _ 2^(£) 

c cz c 


(3.4) 


g(t) - | f ( t) | 2 * h(t) , (3.5) 

K = A^(Xz) 2 [ /d 2 _p |a(j),z)| 2 ] 2 / /d 2 _p | a (_p , z ) | ^ , (3.6) 

<N> " TTF 0 r Q T a 2 \ Z ~ 2 • (3 ‘ 7) 

B° “ / d“£ l fl (£» z ) | 2 B r (_p) / /d 2 £ | a (_p , z ) | 2 (3.8) 


with 

p = (x,y) ■ horizontal coordinate vector on the target surface measured 
from the center of the footprint, 
a(p,z) * complex amplitude cross-section of the laser footprint, 

B r (£) “ power reflection coefficient of the target, 

8° * average power reflectivity of the target within the footprint, 
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£(j>) * surface profile of the target, 

z “ target distance, 

c - velocity of light, 

2 

j f ( t ) | * transmitted laser pulse intensity, 

h(t) - impulse response of the receiver electronics, 

A ” area of the receiver aperture, 

K 

\ - laser wavelength, 

Q ■ total energy transmitted per laser pulse, 

T * one-way intensity transmittance of the atmosphere, 
n “ quantum efficiency of the photodetector, and 
hf ■ energy of one signal photon. 

Throughout this thesis the spatial integrals are assumed to be eva- 
luated over the entire _p plane. The pulse propagation time ( 4>) is a func- 
tion of the surface profile of the target £(j>) within the laser footprint 
which is generally unknown and may be considered to be random. The first 
term in is the nominal round-trip propagation time from the ranging system 
to the center of the target. The second term is the additional delay due to 
the curvature of the laser wavefront, and the last term is due to the range 
spread of the target. <N> is the expected number of detected signal photons 
per received pulse. K is the. signal-to-noise ratio of the speckle which can 
also be regarded as the number of speckle correlation cells averaged by the 
receiver aperture. K is equal to or greater than one and typically is much 
larger than one. The first term in Eq. (3.2) is due to the 6hot-noise, 
while the second term is due to the speckle. 
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3 . 3 Statistics of the Reflected Signals from the Cube Corner Reflector 

(CCR) Arrays 

Most optically rough targets with sufficiently large numbers of 
independent scattering areas will yield the signals statistics that were 
described in Section 3.2. However, there are targets of interest which do 
not fall into this category, and the resultant statistics of the reflected 
signals can be quite different. One important example is the Laser 
Geodyramic Satellite (Lageos) - the first space target devoted exclusively 
to Laser ranging. Lageos is a passive, long-lived target with well defined 
orbit [Fitzmaurice , et a)., 1977]. It is a 60-cm diameter sphere orbiting 


at an altitude of 5900 km. Its surface is covered with 426 cube corner 


reflectors (CCR) to retrodirect any incident optical signal. Because of the 
relatively small number of point reflectors on Lageos, the reflected 
electrical field does not exhibit Gaussian statistics, and the resultant 
speckle is the so-called "partially developed" speckle [Goodman, 1975]. 

The reflected electric field from any CCR array consists of the 
coherent addition of the reflected electric field components from CCR's at 
various T anges within the array. The phase angles associated with these 
field components are in general random. If the range or height difference 
between two or more CCR's is shorter than the transmitted pulse width, the 
reflected pulses will overlap in time at the receiver, resulting in either 
constructive or destructive interference depending on the relative phase 
difference. In this case the detected power at the receiver will vary ran- 
domly from pulse to pulse as the relative phase changes due to changing 
target orientation. This coherent interference can introduce random fluc- 
tuations in the detected pulse shape, and consequently, increase the range 
error [Bufton, et al., 1977]. 
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The mean detected signal and the signal autocovariance for reflections 
from a CCR array can be derived by using the discrete facet model for specu- 
lar reflections, the Fresnel diffraction formula and Campbell's Theorem 
[Papoulis, 1974]. The detailed calculations are shown in Appendix A. The 
results are [Im and Gardner, 1985] 

M 9 

<S(t)> = <N> l B^.z) | f(t - ♦(Pj) I * h(t) , (3.9) 

m=l 


and 


Cg( t ^ , t 2 ) = <N> I J dT I ~ ^))| h ^ c i " T ) h ^ c 2 “ T ) 

m=l -®> 


<N> 2 K _1 l l B 4 (p^.Z.A) J’dTj f(Tj 
m=l n=l 
m*n 


M M 


4>( D ) ) 

—Cl 


•^(Tj - 4'(P n )) h(tj - Tj) / dx 2 f*(x 2 - *(^ n )) 

— 00 


(3.10) 


•f(x 2 - iK^)) h ( 1 2 - t 2 ) 


where 


B (p ,z) = 

L —Cl 


tn 


a< ^»' 


; ) r / 


M 

l 

ra=l 


m 


a(p ,z) |‘ 
— m 1 


(3.11) 
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B 4 ( An’A»’ Z ’ X) 


9 ( ^ 2 

0 o ja(p ,z) j" | a ( o , z ) j .2 . 

ra n 1 An 1 ' — n 1 — m — n 


M M 

l l ° m °n l a( -Bm « z) | 2 l a( Ai’ z) I 2 w 2 ( -Sm*-en' Z ' X) T( An*Ai } 

m=l n= 1 


m*n 


(3.12) 


K = 


M 


l 0 m laCj^.z) | / dt | f ( t ) |‘ 


m=l 


(3.13) 


M M 

I l ° m ° n l a( £d * z) | 2 l a( Ai ’ z) ! 2 ^An’A^*’^ ^An’Ai* 

m=l n=l 
m*n 


<B> - h GC > T l \ * 


(3.14) 


G 


V a |a(p ,z) / Jd^p |a(p,z)|^ 

^ tn ' “Hu 


(3.15) 


y(p, J2 n ) = I ^ dt f(t ‘ * ( 4 b }) f * (t ‘ 


(3.16) 


a 

ra 



W( An’ Ai’ z * X) 


J 1 [2(Xz)' 1 (^y 172 !^ - Pnll 
[( Xz) -i (^V 1/2 La, " p n N 


(3.1/) 


(3.18) 


with 
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p ■ horizontal coordinate vector of the m-th CCR, 

— m 

M * number of CCR's in the array, 

* effective reflecting area of the m-th CCR, 

= power reflectivity of the m-th CCR, and 
m 

G = gain of the CCR array. 

J, is the Bessel function of the first kind of the first order, o is 
1 m 

the effective lidar cross-section when the specular reflection from the m-th 

CCR is redirected toward the receiver [Fitzmaurica et al., 1977]. The 

effective lidar cross-section of the CCR's at optical wavelengths is usually 
6 2 

very large (~10 in ). The effective reflecting area (A ) is a function of 

m 

the orientation of each CCR and is generally smaller than the actual surface 

area. The speckle SHR is dependent on both the optical wavelength and the 

range distance. The coherent interference factor (y(p ,p )) is the time 

— m — n 

integral of the product of two transmitted pulse envelopes that are 

separated by the differential time delay ijKp^) ~ By referring to Eq. 

(3.4), the differential time delay is directly proportional to the height 

variations between two CCR's (£(p ) - £(p )). Therefore, when the two CCR's 

— m — n 

are separated sufficiently in range, YO^.J^) is zero. When they lie on the 

same plane, y(p ,p ) is one. W(m ,o ,z,A) is the aperture weighing factor. 

When the range distance is long so that each individual CCR cannot be 

resolved by the receiving telescope, W is approximately equal to one. 

However, for short distance ranging, individual CCR may be resolved and W 

may be significantly less than one. 

Equations (3.9) through (3.18) are valid for ranging to CCR array at 

any arbitrary distance. However, for most practical applications the range 

3 4 

distances are on the order of 10 -10 km so that the following inequality 
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< < Xz 


V m,n 


(3.19) 


is satisfied. In tnis case, W is approximately 1 and Eqs. (3.12) and (3.13) 
can be simplified to 


B 4 ( 4n»£o» z) =' 


°m °n l a( -£m’ z) ^ l a( Ai‘ z) ^ 


(3.20) 


M M 

l l ° m ° n I z I z 

m=l n=l 
m*n 


and 


K = 


l a m I" /°° dt l f(t > r 

m=l -<*> 


M M 

l l a a |a(p ,z)| Z |a(p ,z) | Z Y(p >P ) 
m=l n=l 
m*n 


— m — n 


(3.12) 


respectively. Therefore, at reasonably long range, the speckle SNR 
for reflection from a CCR array is independent of the optical wavelength, 
but is affected strongly by the coherent interference. In this regard, the 
speckle SNR for reflections from a CCR array is quite different from that 
for reflections from a diffuse target. The speckle SNR for reflection from 
a CCR array and its importance in laser ranging applications will be 
discussed in more detail in Chapter 5. 

3.4 Gaussian Model for the Transmitted Pulse Intensity and Receiver Impulse 
Response 

2 

If the transmitted pulse intensity (|f(t)| ) an 4 receiver impulse 
response (h(t)) are modeled as Gaussian functions of time, <S(t'> and 
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Cg ( c i » 1 2 ) ^ or c ^ e reflected signals from both the diifuse target and CCR 
array can be simplified to give 


<S( t )> = <N> f s (t) , 


(3.22) 


and 


C 1 + C 2 , 


C S (t l ,t 2 ) = <N> G(v ^” V t \ ~ t 2 ) f S ( 2 } 


(3.23) 


?i t + t- 

♦ <N> 2 K G(/2 o tj - t 2 ) f sp ( ^ ) 


where 


2 

)d _p b 2 (_p,z) G(o^, t - ip(_p) ) diffuse target 


f s (t) “ l 


(3.24) 


M 

J B_ ( p ,z) G(o , t - t|/( p )) CCR array, 
2 —m g — m J 

m -1 5 


/ d 2 _p b^(_p,z) G( o / /2 , t - \K_p) ) 


diffuse target 


fgpt 1 7 


(3.25) 


M M 


l l V -2 m’Ai* z) /fi, t 


i|»( p ) + <K P ) 






m=l n=l 
m*n 


) • Y(P ,M_) 

—in — n 


CCR array, 
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Y(P ,p) = exp { 

— m —fi 1 


[<K_P ) " i|>(f> )1‘ 

— m — n 

4a, 2 


(3.26) 


G(a,t) = (/Tn a) ^ exp [-t 2 /2a 2 ] , 


(3.27) 


2 2 2 
a = a + a . 
g h f 


(3.28) 


fg(t) is the normalized waveform of the mean detected signal. 

fgp(t) is the normalized waveform of the speckle induced variance of the 

received signal. Both f g ( t ) and fgp(t) are functions of the target 

geometry. o- and are the rms widths of the transmitted pulse intensity 

( | f | ) and the receiver impulse response (h), respectively. is the rms 

width of the reflected pulse intensity from a nondispersive point target. 

In deriving Eq. (3.23), we have used the fact that che width of the detected 

signal is large compared to both o^ and a^. We can see from Eq. (3.23) that 

shot noise has a correlation length of about 2^1 0^ , and speckle has a 

correlation length of about l/l o^. By referring to Eqs. (3.4) and (3.26), 

y(p ,p ) is significant when the following condition is satisfied: 

- m — n 


>Kp ) - V(p ) 

— m — n 


2 2 
P - P , 

— — + - U(p ) - 5(p )] < 2/2 a, . (3.29) 

cz c — m — n t 


Since the optical beam curvature effect (p /cz) is very 3mall, we have 

m 

equivalently 


[£(p ) - £(pJl < 2 CO. 

—to — n 


(3.30) 
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Inequality (3.30) simply states that only the OCR's that are separated 
in range by less than the width of the transmitted pulse will interfere with 
each other. Therefore, short laser pulses must be used in order to minimize 
speckle. In this case, interference only occurs when two or more OCR's hap- 
pen to lie on the same horizontal plane. Equations (3.22) and (3.23) will 
be used in subsequent chapters to analyze the receiver's timing performance. 
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4. COMPARISON OF THE ML ARRIVAL TIME ESTIMATORS 
IN THE PRESENCE OF BACKGROUND NOISE 


4 . 1 Introduction 

In addition to the effects of atmospheric refraction and turbulence, 
the performance of laser ranging systems is also limited by the presence of 
noise. The problem of estimating the arrival (propagation,) time of laser 
pulses in the presence of both shot noise and background noi6e has been 
investigated by Bar-David [1969, 1 973 J . Under the assumptions that the 
envelope of the transmitted pulse is piecewise differentiable and that the 
speckle noise is negligible, Bar-David derived the Maximum-Likelihood (ML) 
estimator, which consisted of correlating the received signal with the 
logarithm of the mean received signal and searching for the peak of the 
correlation function. This estimator is optimal if the target surface is 
assumed to be optically flat or if the range spread of the target is smaller 
than the laser pulse width. However, for the extended diffuse target where 
the range spread is larger than the transmitted pulse width, the width of 
the received pulse will be longer than the correlation length of the speckle 
induced fluctuation. As a consequence, the time resolved speckle will cause 
random small scale fluctuations within the received pulse and introduce 
additional errors in the estimate of the arrival time. In this case, 
Bar-David's estimator will no longer be optimal. In a recert paper, Tsai 
and Gardner [1985J studied the ML estimator of the pulse arrival times in 
the presence of both shot noise and time-resolved speckle. Their results 
show that, under the background-noise-free condition, the accuracy of the ML 
estimator improves directly with the width of the receiver observation 
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interval when the received signal is large. Thus, their estimator appears 
to perform better than Bar-David's estimator at high signal levels over the 
optically rough targets. However, as the width of the observation interval 
increases, the background noise which enters the receiving telescope becomes 
significant and may affect the accuracy improvement of this new ML estima- 
tor. The following analysis investigates the performances of both ML esti- 
mators on the arrival time estimation in the presence of background noise, 
and the results obtained should be more suitable for realistic ranging 
appl ications . 

4 . 2 Statistics of Signal and Noise 

In a direct-detection laser ranging system geometry, the reflected 
laser pulses are detected by a receiving telescope and the photodetector. 

The detector output is then filtered and sampled periodically at an interval 
T equal to the inverse of the filter bandwidth. Thus the sampled values are 
proportional to the number of photons detected in consecutive time bins of 
width T, and the receiver observations can be modeled as a photocount vector 
k, 


k * (kj ,k 2 .... ,k n ) 


(4.1) 


where k. is the photocount within the i-th time bin of width T. k. in 
l l 

general consists of both the signal photocount (k<. ) and the background 

i 

noise photocount (k^ ). That is, 

i 
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k i “ k S * k b * i_1 n - (4.2) 

i i 


The statistics of kg are related to the signal energy W r received 

i 1 

during the i-th time bin. For reflection from diffuse targets and in the 

presence of fully developed speckle, it has been shown that W is 

u • 

1 

approximately gamma distributed [Goodman, 1973 and 1973], obeying the proba- 
bility density function 


P (Wg ) - (r(M.)Wg f 1 [M.Wg / <W g >] l exp[-M i W s /<W g >] , (4.3) 

i i i i i i 


where <W C > is the mean signal energy in the i-th time bin. M. is the 
1 

speckle s ignal-to-noise ratio in the i-th time bin, which is defined as the 
square of the mean divided by the variance 


M. = <W_ > 2 / var(W ). (4.4) 

1 b • b • 

1 I 


The mean and variance of Wg have been derived by Gardner [1982] for 

i 

reflections from rough diffuse targets. They are given by 



> 


<v 


/ dt / d 2 _p b (_p,z) 
(i-l)T 


f(t - \tfp)) 


(4.5) 


and 
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var(W ) = <Q-> 2 K _1 / dt / dt, / d% b, ( p, r) | f ( t, -*( p) ) 

i (i-l)T (i-l)T 4 “ 1 “ 


iT 


iT 


(4.6) 


| f ( 1 2 ~ 4>(_P) ) 


where <Qg> is the expected received signal energy per pulse. The remaining 
quantities were defined in Chapter 3.2. Without loss of generality, we have 
assumed the bandwidth of the photodetector to be much larger than the band- 
width of the transmitted signals in both Eqs . (4.5) and (4.6). 

For a Gaussian-shaped transmitted pulse intensity with RMS width o^, 
Eqs. (4.5) and (4.6) can be simplified so that 


<Wg > = <Q S > F s (iT)T 
i 


(4.7) 


and 


var (Wg ) 3 <Qg> 2 K _1 F sp (iT)T , 
i 


(4.8) 


where 


Fg(t) - / d 2 _P b 2 (_p,z) G(o f , t - *(p)) (4.9) 

and 

F g p(t) " / d2 -2 C(o f //2, t - <|»(£)). 


(4.10) 
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F^(t) is the normalized waveform of the mean received signal prior to 

sampling and Fgp(t) is the normalized waveform of the speckle induced 

variance of the received signal prior to sampling. They both are functions 

of the target geometry. In deriving Eqs. (4.7) and (4.8), we have made the 

assumption that the width of the time bin T is long compared with the width 

of the transmitted pulse, but it is short compared to the widths of Fg and 

Fgp. For the photodetector with finite bandwidth, o^. must be replaced by 
1/2 

(o, + o.) in both Eqs. (4.9) and (4.10) where a is the R>iS width of the 
t d d 

impulse response of the photodetector. 

By using Eqs. (4.7) and (4.8), NT can be expressed as 

M. = K Fg(iT)T / F Sp (iT) . (4.11) 

may be thought of as the number of speckle correlation cells within the 
receiver aperature during the i-th time bin. For targets where the range 
spread is large compared to the transmitted pulse width, M • will change with 
time as different regions of the target are illuminated by the propagating 
pulse. At the receiver, the speckle correlation area is inversely related 
to the area of illumination on the target. If only one point is illuminated 
on the target, the speckle correlation area is infinite and is one. When 
a large number of points are simultaneously illuminated, the speckle will be 
fully developed and will begin to increase. When the target is comple- 
tely illuminated by the laser footprint, M . will be equal to K. Typically, 
will be approximately one near the leading and trailing edges of the 
received pulse and larger than one near the center of the pulse. For 
targets with large range spread, typically only a part of the target area is 
illuminated within one time bin and M. will be smaller than K. 
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Since W Q is gamma distributed with parameter M. , the distribution of 

t> . 1 

1 

kg can be shown to have a probability density function given by [Goodman, 
i 

1963] 


P (k S. ) 

l 


r(k s. * M i> 

1 

r(k s ♦ i) r(M i ) 

i 


M. 


1 + 


<k S. > 

l 


-k 


S. 

L 


1 ♦ 


< k s.> 

1 


M. 

l 


-M. 

l 


(4.12) 


where 


<k s.> ' rr <u s > ■ <N> F s liI) T 

l l 


(4.13) 


with <N> being defined in Chapter 3.2 as the number of expected photons per 
received pulse. If both k c and M. at each time bin are integers, the function 

b • X 

1 

given by Eq . (4.12) is negative binomial. Therefore, the variance of kg is 

i 

var ( kg ) - <kg > ♦ <k g > 2 /M. . (4.14) 

i i i 

For most optical receivers, the target-reflected laser pulses are usually 
much broader than the transmitted pulses. Therefore, the width of time bin, T, 
is larger than the transmitted pulse width. In this case, the signal 
photocounts from different time bins are statistically independent. The joint 
density function of _kg is therefore 
n 

p(Jig) * IT p(kg ) (4.15) 

i*l i 

and the cross covariance is 

f var ( k_ ) 6 . . i * j 

I S i l * J 

cov(kg ,k ) » \ . (4.16) 

i i I 


Background noise is mainly due to the undesirable background radiation 
entering the receiving telescope. External noise can arise from any background 


35 


sources including the sun, the moon, and the stars. Often a more important 
source of noise may be sunlight reflected directly off the target itself. All 
such backgrounds are thermal in origin. Mandel ( 1959) has sb'. m. that thermal 
background noise produces Poisson distributed noise photocounts when the 
receiver time bin width T is large compared to the background noise coherence 
time. In this case, the background noise probability density function is given 

by k b. 

<k b.> 1 


p(k b ) - ~ ‘ T~ ex Pl-< k b .>) 
10 . 1 


(A . 1 / ) 


l 


where 


<k fe > ° var(k b ) - 


(A. 18) 


and X is the average background noise photocount rate, 
b 

Under the assumptions that the noise photocounts are mutually independent 
at different time bins, and that they are statistically uncorrelated with the 
signal photocount at every time bin, we can derive the second order statistics 
as 


cov(k, ,k, ) 

D . D . 

1 J 


var(k. ) 5. . 

1 » J 


(A. 19) 


i * j 


and 

cov(kg ,k b ) ■ 0 i,j, (A. 20) 

i j 

respectively . 

A . 3 Estimation of Arrival Time 

In laser ranging applications, the arrival time of the reflected pulse, 
denoted by t, is ot interest. In this case, Eq . (A. 12) can be written 
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where arg max [F(x)] denotes the argument x which maximizes F(x). By 
x 

substituting Eq. (4.21) into Eq . (4.22) and making the assumption that the 
received pulse always stays within the observation interval for all values of 
T so that there are no end effects, the estimator given by Eq. (4.22) reduces 
to [Tsai and Gardner, 1985] 

T SPEC = arg max t H(T )l (4.23) 


where 

n , n 

K( t) = l (k + M.(t) - -) *n(k + M. (x)) - l k q in 

. , o. 1 L u. 1 

1=1 i i 1=1 i 

(4.24) 

In deriving Eqs. (4.23) and (4.24), we have made use of the following 
approximat ion 

In T(x) = (x - 1/2) In x - x t 1/2 (In 2w) x » 1 . (4.25) 

In the absence of background noise, Eqs. (4.23) and (4.24) describe an 



optimal estimating scheme. This estimator emphasizes signal counts in the 
tails of the received pulse where the speckle is less severe. Therefore, in 
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the large signal limit, the accuracy of this estimator improves as the width of 
the observation interval (nT) increases. However, when the background noise is 
present, its effects will be significant in the tails of the received pulse 
where the signal counts are low. Consequently, we would expect the improvement 
on the performance of this estimator to be offset somewhat by the effects of 
background noise when the observation interval is large. 

In the presence of background noise, Eq . (4.24) becomes 


n . n 

H(t) = I (k. + M.(t) - rr) An(k. + M-(x)) - Z k- °n 
i=l 1 1 L 11 i“l 


M i ( l) 

1 + " <k s" (x)> J 

1 (4.26) 


where ki is given by Eq . (4.2). 

Notice that by replacing Eq. (4.24) by Eq . (4.26) the estimator has become 
subopt imal . 

A 

Assuming that the maximum value that has been chosen using T gp££ is in the 
correct neighborhood of the actual arrival time, T q , the variance of this 
estimator can be calculated in terms of the first two derivatives of H( t) , 
evaluated at the actual arrival time of the pulse. This method of calculation 
has been used in similar problems [Helstrora, I960]. If H( t) is twice 
differentiable about t q , then H( t) can be expanded by Taylor series to become 

H(t) - H( T ) + (t - T ) H( T ) + ... (4.27) 

o o o 


where the dot ( •) has been used to denote the time derivative. By the 
definition of ' l 5 p EC > 


0 h(t spec ) ~ H( V 4 ( t spec " V H(T o ) 4 


(4.28) 


from which it follows that the error of this estimator is 


I 

i 


t spec " T o 


H <V 


H( V 


(4.29) 


provided H(t q ) * 0. This constraint excludes pulses which have constant slopes 

such as the rectangular and triangular pulses from the analysis. 

If we assume that the shot noise and speckle are not severe so that 
2 

<H(t )> >> var[H(t )], then the bias and variance of this estimator can be 

o o 

written as 


<T - x > = - <H(t )> / <H( T )> 
SPEC o o o 


(4.30) 


var(T SPEC ) = var[H(l o )] 1 <H(T o )y 


(4.31) 


In these equations, the expectations are over the observation vector _k. 

By using Eqs. (4.13) through (4.20), it can be shown that <H(t q )> is zero. 

Therefore, t estimator is unbiased. The variance of H( x ) can be 
Sr EC O 


calculated to give 


<k b Xkg > 
i i 



f < k s>l 

|2 

1 

<k > z 

1 ♦ 1 


s. 

M. 
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1 1 J 



(4.32) 


<k S > n <k b ><k S > 

var [ H( i ) ] - l » <k l \ k ' >^ 2 ' (4>32) 

0 i=l ^ k S. i=l I K S/ 

^ 1 + -TT- <k s.> j 1 

1 l 1 J 1 1 1 J 

The first term in Eq. (4.32) is the variances due to the combined effects 

of the shot-noise and speckle in the background-noise-free situation; the 

second term is the additional error due to the presence of background noise. 

<H( t )> can be calculated to give 
o 


<H( x )> = - ) 

° i=l 


<k > 

l 



<k s.>] 

<k S.> 

l 

1 + 

M i j 


(4.33) 
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Notice that <H(t q )> Joes not depend on the background noise counc (<K^ >). 

i 

Finally by substituting these results into Eq . (4.31), we have for the variance 


of the T 


SPEC 


estimator 


Var<T SPEC ) " °T 


SPEC. 


n 

l 

i=l 


<k s.> 
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<k s.> 

1 


1 + 
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M. 
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The first term in Eq. (4.34) is the MS timing error in the background- 
noise-free situation. The second term arises due to the presence of background 
noise . 


As a comparison, the ML estimator derived by Bar-David [1969] considered 
both the shot noise and background noise. Therefore, it is optimal in 
situations where speckle noise is negligible. Mathematically, it can be 
written as 

n 

t^hqt 3 ar 8 max i k i ^n<k i (T)> . (4.35) 

T i= 1 

By employing the same procedure used above, we find it to be unbiased, and 
the variance in the presence of speckle is given by 
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When the speckle noise is small compared to shot noise (i.e., when 
K >> <N>), and when there is no background noise (i.e., ^ = 0), both 
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4 . 4 Timing Performance for Diffuse Flat Target 

To illustrate quantitatively the effects of background noise, we shall 
look at some specific target configurations and evaluate their corresponding 
timing variances. In the first example we consider pulse reflections from an 
infinite flat diffuse target. We also assume that the target surface is normal 
to the ray of laser incidence (Fig. 4.1 with $ = 0). When the laser beam is 
incident normally onto the target surface, the broadening of the received pulse 
is due to the wavefront curvature of the laser beam. 

To simplify the mathematics, we assume that the power reflectivity is 
uniform over the entire footprint, and that the laser footprint has a Gaussian 


cross-section 

| a (_P , z ) | ^ * Q(2nz i tan 2 9^,) * exp(-p~/2z^ tan 2 6^) (4.38) 

where Q is the transmitted laser pulse energy, and 8 is the RMS laser beam 
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divergence angle. In this case, the specxle SNR for the entire observation 
interval (K) can be calculated using Eq. (3.6) to give 


K = tiA r (2 tan 0 T / X)‘ 


(4.39) 


For normal incidence, the mean received signal and the waveform of the 
speckle variance at each time bin, Fg and Fg p , can be calculated to give 
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respectively. Equations (4.40) and (4.41) are plotted in Figs. 4.2 and 
4.3, respectively, for three different laser beam divergence angles. We can 
see that the leading edges of F g and Fg p are sharp and undistorted . However, 
the trailing edges are broadened by the beam curvature effect. 

The speckle SNR at each time bin (Mi) is given by 
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Fig. 4.2. Mean received waveform for normal incidence on an infinite 

flat diffuse target, (z = 500 kra, cCf = 0.5 cm, Ap = 100 ctn^ , 
X = 1 . 064 ym) . 


00 u 



SPECKLE VARIANCE 



TIME 


U-2Z/C) 


cr. 


f 


Fig. 4.3. Speckle variance for normal incidence on an infinite flat 
diffuse target (z = 300 km, cOf * 0.5 cm, Ar * 100 cm^ , 

X = 1 . 064 urn) . 
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Fig. 4.4. Speckle SNR for normal incidence on an infinite flat diffuse 

target, (z 3 500 km, cCf - 0.5 cm, A R - 100 cm*- , X = 1.064 pm). 
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Fig. 4.4. Speckle SNR for normal incidence on an infinite flat diffuse 

target, (z * 500 km, ca f - 0.5 cn, A R - 100 car , X * 1.064 pm). 
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Fig. 4.7. Comparison of the RMS ranging errors using t sp£c and T SHOT 
estimators for normal incidence on an infinite flat diffuse 
target, (z * 500 km, 9- • 0.5 mrad, co^ ■ 0.5 cm, A R - 100 cm 
X * 1 .064 ym) . 
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T 'PEC e8t ^ mator * However, at high signal level where the speckle dominates, 

A A 

the T gp££ estimator performs considerably better than the Tg^^ estimator. 
Consequently, in applications where the signal level is expected to be low, 

A 

T SH0T e8t ^ mator should be used to estimate the pulse arrival times. On 
the other hand, in applications where speckle is expected to be the dominated 

A 

factor, T gpE£ must be used to minimize the timing error. 

For non-normal incidence on an infinite flat target the broadening of the 
received signal is dominated by the range spread of the target. This range 
spread arises because the target in this case is tilted with respect to the 
optical axis (refer to Fig. 4.1). By neglecting the beam curvature effects, 

Fg and Fg p are given approximately by 


Fg(iT) = G(o t , iT - 2z/c) 


(4.44) 


and 


F sp (iT) = G(o t // 2, iT - 2z/c) 


(4.45) 


where 


2 2 2 

^ 4z tan 9 tan 

l°f + c 


L /2 


(4 ’,6) 


$ is the incidence angle of the laser beam with respect to the target surface 
normal. From Eqs. (4.44) and (4.46), both the mean received signal and the 
speckle variance are Gaussian in shape with the RMS width 0 ^. The pulse width 
consists of the transmitted pulse width (o^) and the effects of the range 
spread of the target 2z tan 9^ tan $ [Gardner , 1982 ] . The value of in this 
case is given by 


M i - K T / 2/7 a T . (4.47) 

Again, is constant with time (i). Equations (4.44) and (4.45) are plotted 
in Figs. 4.8 and 4.9 for three different incidence angles. By comparing 
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Fig. 4.8. Mean received waveform for non-normal incidence on an 

infinite flat diffuse target, (z = 500 km, 8-p = 0.1 mrad, 
C0f = 0.5 cm, Ar = 100 cm , X = 1.064 ym) . 
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Fig. 4.9. Speckle variance for non-normal incidence on an infinite flat 
diffuse target (z = 500 km, 9-p = 0.1 mrad, c Of = 0.5 cm, 

Ar * 100 cra^ , X 3 1.064 pm). 
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Fig. 4.8 with Fig. 4.2, we can see that the received signal in the non-normal 
incidence case is much broader in width than that in the normal incidence 
case. Consequently, we would expect the timing error to be larger in the 
tilted target configuration [Tsai and Gardner, 1985]. 

a 

The variance of the T_ ppr estimator is given approximately by 


2 ° T i \ 4 8 2 (<N>,K) 

’ T SPEC 5 TUT gj«N> ,K) * -^T g 2 (<N>>K) 


(4.48) 


where 
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2 0 T (1 + a <N> K 1 e y 7 ) 


(4.49) 


nT/2a T is the ratio between the width of the observation interval and the RMS 
full width of the received signal. Consequently, we would expect this ratio to 
be on the order of three in order for the receiver to sample to entire pulse. 
Because the MS timing error given by Eq. (4.48) depends on g^ and g^ , 
which in turn depend on this ratio, the timing error will depend on the width 
of the observation interval (nT). Because nT also determines the amount of 
background noise allowed to the receiver, so it is a very important factor in 

A 

determining the timing performance of the T estimator. 
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Equation (4.48) is plotted versus <N> in Fig. 4.10 for nT * 14 and 
for three different values of in order to illustrate the effect of 

A 

background noise on the performance of the x^^ estimator. In Fig. 4.11 
Eq. (4.48) is plotted for three different values of nT/'2o^ in order to show 
the significance of the observation interval width at both the low and high 
signal levels. 

In the shot noise limited region, that is , when <N> « K, Eq. (4.48) can 
be simplified to 
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Because a ^ 
x 


According to Eq. (4.51), the timing error by using T S p E £ estimator will 
decrease as the signal strength increases at low signal level. In the absence 
of background noise, the timing performance is only limited by the pulse width 
(or the bandwidth) of the received signal. But when the background noise is 
present, the timing performance will be degraded inv.rsely with both the 
background noise count rate and the observation interval width. These effects 
can be seen from both Figs. 4.10 and 4.11 in the small signal region. 

is straightly increasing with respect to nT, and therefore, the 
observation interval width must be chosen to be as small as possible in order 
to minimize the timing error. It should be noted that the observation interval 
cannot be chosen to be smaller than the pulse length of the received signal 
( ~ 6 o^.) . When nT is smaller than 60 ^., Eq. (4.51) is no longer valid. However, 
it is easy to show that the timing error in this case will increase as the 
observation interval width is reduced. 

In the speckle limited region, that is, when <N> >> K, the MS timing error 
using T S p EC i- s bounded below by 
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Fig. 4.10. RMS ranging error using t spec estimator for non-normal 

incidence on an infinite flat diffuse target, (z =500 km, 
9-p = 0.1 mrad, $ * 0.5°, caf = 0.5 cm, = 100 cm 2 , 

X = 1.064 ym, nT/2ox = 7). 
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Fig. 4.11. Dependence of RMS ranging error on the observation interval 
width by using T_ ppr estimator and for non-normal incidence 
on an infinite flat' diffuse target, (z = 500 km, = 0.1 mrad 
= 0.5°, cOf = 0.5 cm, Ar = 100 cm , X = 1.064 pm). 
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<N> » K 


(4.52) 


When both and nT are small, the first term in Ineq. (4.52) dominates, and the 
timing error will decrease as the observation interval width increases. 

However, as nT is increased, the exponential term will eventually dominate such 
that the second term exceeds the first term. In this case, any further 
increase in nT would result in the increase of timing error. This effect is 
illustrated by Fig. 4.11 at the large signal region. By applying L'Hopital's 
rule, it can be shown that 
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Consequently, it is theoretically possible for the timing error of T gpgQ 
estimator to approach infinity when the observation interval width is 
exceedingly large. In general, an optimal width can be determined by using 
conventional optimization technique. 

In Fig. 4.12, the RMS timing error is plotted versus nT/2o T for 
different values of <N> and for the typical satellite laser ranging parameters. 
According to this figure, the optimal width of observation interval is between 
4 and 6, this corresponds to the full width of a Gaussian pulse at the 5% 
point. Consequently, nT must be chosen to be not more than the pulsewidth of 
the received signal in order to obtain the optimal performance. Fortunately, 
the rate of increase of the timing error with respect to the observation 
interval width is relatively slow, and therefore, we would not expect 


significant reduction on the timing performance of T gp£ C when the chosen 


observation interval width is a few times larger than the optimal width. 


or 



0 10 20 30 40 5 

5 15 25 35 45 

NORMALIZED WIDTH OF OBSERVATION INTERVAL (- 


Fig. 4.12. RMS ranging error versus width of observation interval using 
tspec estimator for non-normal incidence on an infinite flat 
diffuse target, (z * 500 km, 0-j. = 0.1 mrad, $ = 0.5°, 
cot = 0.5 cm, A R * 100 cm^, X = 1.064 ym) . 
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The timing variance of the estimator can be calculated using Eqs 

(4.36), (4.44) and (4.45) to give 
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Equation (4.54) is plotted versus <N> in Fig. 4.13 for different values 
of It is also plotted in Fig. 4.14 for different values of nT/2^; 

In the shot noise limited region, Eq. (4.54) reduces to 


4/ " *b 


SH0T ' <N> 2 


<N> « K 


(4.57) 


Therefore, in the low signal limit, the timing error for the Tg H ^ T estimator 
will increase directly with the background noise photocount rate. However, Eq. 
(4.57) does not depend on nT, and therefore, will not change even if the 
observation interval width is increased. Tnese effects are clearly shown by 
Figs. 4.13 and 4.14 at the small signal level. Also, this MS error has a 

-2 -i 

<N> dependence instead of the <N> dependence as in the background-noise- 


free situation. 
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Fig. 4.13. RMS ranging error using TgnoT estimator for non-normal incidence 
on an infinite flat diffuse target (z_ a 500 km, * 0.1 mrad, 

$ * 5*, cOf ■ 0.5 cm, Ar * 100 cm^ , \ * 1.064 urn, nT/2c>p “ 7). 
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Fig. 4.14. Dependence of RllS ranging error on the observation interval 
width by uaing estimator and for non-normal incidence 

on an infinite flat diffuse target, (z * 500 km, 0y 3 0.1 mrad, 
$ - 0.5" 1 , cof 3 0.5 cm, Ar ■ 100 cuT , X = 1.064 pm). 
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In Che speckle limiting region, Eq. (4.54) reduces Co 


SHOT 


2K 


<N> « K . (4.58) 


From Eq. (4.58), Che Ciming performance of T g. rJ0T does noC depend on eicher Che 
background noise or Che observacion inCerval width. In Chis region, however, 
Che speckle effecCs become significanC. Consequently, Che performance of Chis 
escimator will saCurate at high signal level. 

To compare Che performance of Chese Cwo estimators on the tilted diffuse 
flat Carget, both Eqs. (4.48) and (4.54) are plotted in Fig. 4.15 for 
b ■ 3xl0^/sec. The results show that within normal background noise level, 

A A 

the performance of the T gjjQj estimator is better than that of the T gpg£ esti- 
mator at low signal level. However, at high signal level, the T gpgc estimator 
will perform better than the T gnoT est i mator * 

4 . 5 Discussion 

In laser ranging applications, snot noise, speckle and background noise 
introduce ranging errors, thus, limit the ranging (timing) accuracy of the 
receiver. In this chapter, two Maximum- 1 ike 1 ihood arrival time estimators were 
investigated in order to understand their performances under the influence of 
the background noise. Bar-David's t shqt estimator takes into account 
both shot noise and background noise, and therefore, its performance is 
relatively less sensitive to the addition of background noise. However, at 
high signal level, its performance is limited by the time resolved speckle. On 
the other hand, Tsai and Gardner's T spi - C estimator considers both shot noise 

A 

and speckle. Consequently, it performs better than at high signal level 

where the speckle is the dominated error source, but it is subject in a greater 
extent to the addition of background noise. Also, this estimator is dependent 
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Fig. 4.15. Comparison of the RMS ranging errors using igpE^ and 

estimators for non-normal incidence on an infinite flat 
diffuse target, (z * 500 km, 9-j. - 0.1 mrad, $ - 0.5°, 
cOf * 0.5 cm, A r * 100 cm^ , X • 1.064 um) . 
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on Che width of the observation interval, especially at low signal level. At 
high signal level, its performance will be degraded slightly when the width of 
the observation interval is large compared to the pulse width of the received 
signal . 

Background noise effect appears to be most sigificant at low signal level in 
the estimation of the arrival time. At high signal level, however, its effect 

A 

is negligible. In normal background noise level, Tg H Q T estimator will perform 

A A 

better than T gp EC estimator at low signal level, but T gp EC estimator will 
perform better at high signal level. 
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5. TWO-COLOR CROSS-CORRELATION ALGORITHM 


5 . 1 Introduction 

In satellite laser ranging applications where it is impossible or too 
expensive to obtain accurate meteorological information around the ranging 
sites, or in applications where optical pathlengths at two frequencies are 
necessary in order to obtain the appropriate parameter information, such as 
the surface pressure, it is desirable to measure the difference in the 
round-trip propagation time between pulses that were transmitted simulta- 
neously at two optical frequencies. The Maximum-Likelihood estimators 
described in Chapter 4 can be used to determine the differential arrival 
time of the detected pulses in two-color ranging. Unfortunately, implemen- 
tation of the ML estimators requires a prior knowledge of the mean detected 
pulse shape, and therefore, can only be used in applications where the 
target profiles and orientations are known. In most cases of interest, 
however, the mean pulse shapes will not be known and suboptiraal estimators 
must be used. 

The centroid estimator is easy to implement and is optimum for Gaussian 
pulses when speckle is negligible. Unfortunately, most of the reflected 
pulses from realistic targets are non-Gaussian and asymmetric, which reduces 
the effectiveness of the centroid estimator. Iyer [1976] has proposed a 
leading-edge threshold detection technique which registers the arrival time 
the instant that the output of the optical detection system exceeds a 
preset threshold. Receivers of this type are frequently used in laser 
ranging applications. However, this algorithm introduces a signal dependent 
bias which must be compensated in high precision ranging applications. 
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In this chapter a two-color cross correlation technique is being 
studied; the results will be applied to the pulse reflections from ground 
diffuse targets, earth-orbiting CCR arrays, and ocean. 

5.2 Cross-Correlation Algorithm 

In a two-color ranging system, the laser pulses of t»*o different 
wavelengths are transmitted simultaneously. If the two beams are aligned so 
that their footprints overlap, the mean reflected pulse shapes will be 
almost identical. An effective technique for estimating the differential 
arrival time that takes advantage of this pulse similarity is calculation of 
the peak of the cross-correlation function of the two detected signals. 

This correlation estimator can be written mathematically as [Gardner et al., 
1983] 



arg max 
T 


/ dt S^t) S 2 (t + T) 


arg max [R 12 ( » (5.1) 

T 


where and S 2 are the detected signals at wavelengths X, and X 2 , 
respectively, Rj 2 (t) is the cross-correlation function of these two signals, 
T is the differential arrival time, and is its estimate using the 

correlation algorithm. 

To analyze the performance of the correlation estimator, we must 

linearize R,_ around the actual differential arrival time ( t ) by using a 
12 o 

Taylor series expansion. This method depends upon the underlying assumption 
that Rj 2 is differentiable. Therefore, we shall exclude from the analysis 
detected signals which would give rise to a non-dif ferentiable R^ 2 such as 
reccangular pulses. With the assumption that the shot noise and speckle are 


not severe so that 
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var lR 12 ( T 0 )] < < E 2 [ r 12 ( t 0 )) 

the bias (E(Ax )) and the variance 
uUK 

can be calculated to give 
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Eq. (5.4) gives the lower bound of the correlation estimator. The actual 
timing variance is expected to be higher depending upon the actual SNR. 

When the SNR is high, however, the actual timing variance will approach this 
bound asymptotically. 

By noting that both the speckle and the shot noise induced fluctuations 
are uncorrelated at the two frequencies, we can express the bias and the 
variance of the correlation estimator in terms of the mean and autocorrela- 
tion functions of S^ and S 2 . That is, 
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where 

R s (t 1 ,t 2 ) - Cg ( c i * 1 2 ) + <S i (t 1 )XS i (t 2 )> , i-1,2. (5.7) 

l i 

The mean (Eq. (3.22)) and autocovariance function (Eq. (3.23)) at the two 
frequencies are related by 

<S i (t)> = <N i > f s (t - (i-l)T o ) i-1,2 (5.8) 

and 


1 1 + t 2 
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(5.9) 

, t. + t, 

+ <N.> Z K. i G(/2 a , t l - t 2 ) f sp ( ^ - (i-l)T Q ) , i-1,2, 

respective! y. By substituting Eq. (5.8) into Eq. (5.5), it is easy to show 
that the correlation estimator is unbiased. Since this estimator cross 
correlates the two noise corrupted signals, the peak of the cross correla- 
tion function will shift randomly and introduce a timing error. In order to 
analyze the characteristics of this timing error, Eq. (5.6) is expanded 
using Eq. (5.7) to give 
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1 2 -" 1 2 


VV T o 


{ /" dt <Sj(t)XS 2 (t + t q )> } 


The first term on the right side of Eq. (5.10) equals the square of the 
bias, so that it is zero. The second and the third terms are the first- 
order errors involving the cross-correlations of the mean of and the 
fluctuating portion of S^, and vice versa. The last term is the second- 
order error involving the cross-correlations of the fluctuating portions of 
both Sj and S 2 . Notice that the first-order errors will be the total error 
for an estimator that correlates the detected signals with their 
corresponding means. From this point of view, the second-order error is an 
additional error incurred when the mean pulse shapes are not known, which is 
the case in most applications. 
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Because of the random nature of most of the target surface profiles, 
fg and fgp are generally very difficult to evaluate, and numerical com- 
putation is usually required to determine the timing error given by Eq. 
(5.10). It is more convenient to express the timing error in the frequency 
domain. By substituting Eqs . (5.8) and (5.9) into Eq. (5.10) and Fourier 
transforming the various terms, the mean-square (MS) timing error can be 
expressed in terms of the bandwidth of the detected signal. The result is 
[Im and Gardner, 1985] 


a 2 = + IL + J_ + J_j 

T C0R <N 1 > <N 2 > 2B^ [ K 1 K 2 j 2B^ 



where 
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(5.12) 


(5.13) 


(5.14) 


(5.15) 


(5.16) 


(5.17) 


(5.18) 


(5.19) 


(5.20) 
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Bg is the RMS bandwidth of the mean detected signal, Bg p is the RMS hand- 
width of the speckle variance waveform, is the RMS 'cross-bandwidth' bet- 
ween the mean detected signal and speckle variance. Since fg and 
fgp involve the convolution of the transmitted pulse and the receiver 
impulse response, the RMS bandwidths Bg, Bg p and B^ are always larger than 
0 ^ . The ar's are dimensionless factors which typically vary between 1 and 
2 depending on the surface profile. 

In Eq. (5.11), the terms involving <N^> * are the first-order shot- 
noise contributions to the timing error due to the correlation between the 
mean signal at one frequency and the shot-noise induced fluctuating com- 
ponent of the detected signal at the other frequency, and vice versa. The 
terms involving * are the first-order speckle-noise contributions to the 
timing error due to the correlation between the mean signal at one frequency 
and the speckle induced fluctuating component of the detected signal at the 
other frequency, and vice versa. The term involving (<Nj><N 2 >) * is the 
second-order shot-noise contribution due to the correlation between the 
shot-noise induced fluctuating components of the detected signals at the two 
frequencies. The term involving (K^Kj) ^ is the second-order speckle-noise 
contribution due to the correlation between the speckle induced fluctuating 
components of the detected signals at the two frequencies. The last two 
terras involving (<Nj>K 2 ) 1 and (<N 2 >Kj) 1 are the second-order errors due to 
the correlation between the shot-noise induced fluctuating component of the 
detected signal at one frequency and the speckle induced fluctuating com- 
ponent of the detected signal at the other frequency, and vice versa. 

The first-order errors depend on the bandwidth of the mean detected 
signal, which is related to the target surface profile and to the beam 
curvature. The second-order errors depend on the difference between the 
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receiver bandwidth and the bandwidths that are related to the detected 
signal. It should be emphasized that each of the differences in the second- 
order terms are strictly positive. When the receiver bandwidth is very 
large, the timing error is dominated by the second-order shot-noise effect. 
On the other hand, when the receiver bandwidth is very small, the bandwidth 
of the mean detected signal wil? also be very small. Consequently, the 
first-order errors will dominate the overall timing error. In these two 
extreme cases, however, the overall timing error is likely to be very large. 
Fortunately, the first-order errors are less sensitive to the change in 
receiver bandwidth than the second-order errors, and the second-order errors 
can be reduced to negligible levels by reducing the receiver bandwidth 
before the first-order errors diverge significantly. In general, there is 
an optimum receiver bandwidth which minimizes the overall timing error for a 
given set of received pulses. 

The timing error of the correlation receiver is inversely proportional 
to the signal bandwidth. This is expected since high bandwidth signals will 
contain fine scale structure which improve the performance of the correla- 
tion timing algorithm. At low signal levels, the timing performance is 
dominated by either the first-order or the second-order shot-noise effect 
depending upon the difference in the bandwidths of the receiver and the mean 
detected signals. In either case, the timing error will decrease as the 
received signal strength increases. At high signal levels, however, the 
speckle becomes significant, and it places the fundamental limit on the 
timing accuracy of the correlation receiver. 
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5 . 3 Pulse Reflections from Infinite Diffuse Flat Targets 

To illustrate quantitatively the timing performance of the correlation 
receiver and its dependence on the target profile, we shall look at specific 
target configurations. In the first example we consider pulse reflections 
from an infinite flat diffuse target. This is a good model for reflections 
from ground targets. When the laser beam is incident normally on a flat 
surface, the reflected pulses are broadened by the wavefront curvature of 
the laser beam. However, this effect is usually small for typical beam 
divergences. When the laser beam is incident at an angle $ with respect to 
the surface normal (refer to Fig. 4.1), a range spread will be introduced 
by the diffuse target. The amount of range spread for a given target 
depends on $ and the illuminated surface area. For an infinite flat target 
at a range z from the ranging system, and for laser beam divergence angle 9^, 
the RMS range spread is approximately 2 z tan 9^ tan <J> [Gardner, 1982]. The 
range spread induced broadening of a received pulse can significantly reduce 
the signal bandwidth, and, as indicated in Section 5.2, will degrade the 
timing performance of the correlation receiver. 

To simplify the mathematics, we assume that the power reflectivity is 
uniform over the entire footprint, and that the laser footprint has a 
Gussian cross-section 

|a(_p,z) | 2 = Q(2n z 2 tan 2 9 T > 1 exp(-p 2 /2z 2 tan 2 9 T ) (5.21) 

where Q is the transmitted laser pulse energy. fg and f sp can now be 


calculated by using Eqs . (3.3), (3.4), (3.24), (3.25) for the diffuse 
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target. By neglecting the curvature effects, we obtain 
f s < t) - G(o s , t - 2z/c) 

and 


fgp( t) 


G(o s //2, t - 2z/c) 


with 


(Oj + 


2 

°h 


a])* 


<o2 g 


o 2 ) 112 


( 5 . 22 ) 


( 5 . 23 ) 


( 5 . 24 ) 


1 2 z tan 0 p tan $ /c . (5.25) 

Og is the width of the reflected pulse from the flat diffuse target for the 
case of non-normal incidence. o is the range spread induced pulse 
broadening. f g and f gp are plotted in Figs. 5.1 and 5.2 respectively for 
the nominal ranging system parameters and three different incidence angles. 
We can see that both f g and fg p are significantly broadened even for a very 
small incidence angle. The number of speckle correlation cells can be 
calculated by using Eqs. (3.6) and (3.21): 

K . “ir A„(2 tan 0 /X.) 2 . (5.26) 

1 K 11 

Typically, is of the order of a few hundred or more. The MS timing 
error can be calculated by substituting the Fourier tranfonns of Eqs. (5.22) 
and (5.23) into Eqs. (5.11) through (5.20). 


I 


The results are 
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a f 


Fig. 5.1. Mean detected signal waveform (f s (t)) for non-normal Incidence on an 
Infinite flat diffuse target, (z • 500 kn, 0j ■ 0.1 mrad, 

°f ■ 16.7 psec, " 100 psec). 
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Fig. 5.2. Speckle variance (fsp(O) for non-normal incidence on an infinite flat 
diffuse target, (z ■ 500 km, S-p ■ 0.1 mrad, o f - 16.7 psec, 

- 100 psec) 



78 


°t 2 • — ( ?rr* <rr)°s * ) « 

COR 3/T <N 1 > <N 2 > S K 1 K 2 ' ' 


5 , , 2. 2. 5., 2,2. 

t g s (1 ~ q h /2g S ) , g s (1 ~ V°S ) 


/2 <N 1 ><N 2 >0^ 


K,K 0 o 

1 2 g 


(5.27) 


5 °8 * V 


°s (1 - 


• - ( ♦ 1 1 

n 1<N 1 >K 1 <n 2 >k i 


3«; 


<o 2 ♦ „y /2 

g n 


From Eq. (5.27), we can see that the first-order timing errors of the 

correlation receiver increase directly with a g , which in turn is directly 

proportional to the range spread of the target o^, the transmitted pulse width 

Op and the width of the receiver impulse response 0^. The second-order 

errors are i • 'erseiy related to o, . Since a„ > o > a, , the second-order 
J h S g h 

errors will iways exist in this receiver. The timing performance for this 

configuration is plotted against the expected signal strength 

<N^> * <N 2 > = <N> in Fig. 5.3 for the same ranging parameters used in Figs. 

5.1 and 5.2. The plotted result shows that the timing error can be quite 
large due to the range spread effect. The curves saturate at high signal 
levels because of the speckle. 

At low signal levels, the timing error of the correlation receiver is 
dominated by the shot-noise contributions, and Eq . (5.27) can be written 


approximately as 


RMS TIMING ERROR IN PSEC 



NUMBER OF PHOTOELECTRONS (<N>) 


RMS differencial timing error versus received signal strength 
for pulse reflections from an infinite flat diffuse target using 
the correlation estimator (z = 500 km, ftp = 0.1 mrad, 

Of 3 16.7 psec, 0^ = 100 psec, = 0.355 urn, X 2 = 1.064 um, 

Ar 3 0.1 m^ ) . 


Fig. 5.3. 
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V 2B » - B s> 

4/2 <N 1 XN 2 > Lg 


<N . > « K. (5.28) 
1 1 


where 

B s - </2 o s r‘ - [2 <o 2 £ . c 2 h ♦ a 2 r ) f 1/2 
B r = (/2 a ,)' 1 . 


B is the RMS bandwidth of the mean detected signal, and B is the RMS 
receiver bandwidth In this region, if the receiver bandwi^h is large, the 
second-order shot-noise effect is the most significant error source, and the 
timing error of the correlation receiver can be very large. On the other 
hand, when the receiver bandwidth is small, the first-order shot-noise 
effect dominates, and the timing error of the correlation reciever in this 
case is likely to be smaller. In Fig. 5.3, we chose the parameters so 
that B » B_ . Therefore, the RMS timing error curves decrease according to 

K u 

<N> _1 when <N> « K. . 

1 

At high signal levels, the speckle noise dominates, and the timing 
error of che correlation receiver can be written approximately as 


n ( p ^ _ „2 % 

2 1 r \ . 1 1 1 . PT vb PT V 

°r 3 — ir r i ~~o 5 — 

COR /2 1 ^2 Bg 2 K i K 2 B S 


<N { > » K £ (5.29) 


where 
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b pt * (/2 ° g ) ' i " (2(o h + °f )] V2 


(5.30) 


Bp^. is the RMS bandwidth of the point target response of the system. In 

this region, the shot-noise effect is negligible. This behavior is 

illustrated by the saturated (flat) regions of the three curves shown in 

Fig. 5.3. As we can see, the speckle noise places a fundamental limit on 

the accuracy of the differential arrival time estimation. Unlike the 

second-order shot-noise error, the second-order speckle-noise error will 

saturate at a finite level even if B R approaches infinity ( -*-0). This is 

due to the fact that the correlation length of the speckle-induced flue- 

2 2 1/2 

tuation in the received signal is 2/2(a + o^) ' (refer to Eq. (3.23)), so 

even if o^ is very small, there is still a small amount of correlation left 

on the speckle-induced fluctuating components of the received signals, and 

this prevents the timing error from diverging. 

The effect of receiver bandwidth on the timing performance of the 

correlation receiver can be illustrated by plotting in Fig. 5.4 the timing 

error versus the normalized receiver bandwidth (B R /Bg) for three different 

5 « 2 2 - 1 / 2 . 

values of <N> and for = 9 K 2 = 10 ■ = [2(o^- + ) ] is the RMS 

bandwidth of the received signal prior to detection and therefore does not 
depend upon the receiver bandwidth. The upper curve (<N> = 500) corresponds 

A 

to the case of shot-noise-limited differentiil timing error. When B R >> B<, , 
the second-order shot-noise effect dominates, and the overall timing error 
increases with increasing receiver bandwidth. Consequently, the timing 
performance can be improved by reducing the receiver bandwidth. On the 
other hand, when B << B„, the first-order effect dominates, and the timing 
error decreases as the receiver bandwidth increases. We can see that the 


RMS TIMING ERROR IN PSEC 


I0 5 b 


K, = 9K 2 = 10' 


< N > = 500 


I I I I III! I I I I I I III I 


I0' 1 10° I0 1 


NORMALIZED RECEIVER BANDWIDTH (B„/B s ) 


Fig. 5.4. RMS differencial timing error versus normalized receiver bandwidth 
(Bg/Bg) for pulse reflections from an infinite flat diffuse target 
using the correlation estimator. Br is the receiver bandwidth, and 
Bg is the signal bandwidth prior to detection. (z = 500 km, 

9-p = 0.1 rarad , Gf = 16.7 psec, $ = 0.5°, Xj = 0.355 urn, 

X 2 = 1.064 pm, Ar = 0.1 m^) 
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minimum timing error occurs somewhere between these two extremes. The 
middle curve (<N> * 2x10^) corresponds to the case where the shot-noise and 
speckle contributions to the timing error are almost equal. In this case 
the dependence of the timing error on is almost the same as for the upper 
curve. The bottom curve (<N> = 5x10^) corresponds to the timing error in 
the speckle-limited region. This curve has almost the same basic charac- 
teristics as the other two curves except when the receiver bandwidth is very 

A 

large. When >> Bg, the timing error saturates because the second-order 
speckle-noise error is much less sensitive to the change of B^ for reasons 
discussed earlier. 

To illustrate the effects of speckle in the differential arrival time 

estimation, the RMS timing error is plotted versus the signal strength <N> 

for different s>e«.s of ' s in Fig. 5.5. The flat regions in this figure 

correspond to the speckle-limited regions. We can see that the timing error 

increases as the speckle SNR decreases. When the K^'s approach infinity, 

the ranging system operates in the shot-noise-limited region and so the 

timing accuracy is always improved as the received signals become stronger. 

At low signal levels, the second-order shot noise is the most significant 

error source so the curve has a <N> ^ dependence. At high signal levels, 

the first-order shot noise is the most significant error source, so the curve 
- 1/2 

has a <N> dependence. However, for all finite values of K^'s, the 
timing accuracy will saturate when <N> = K^. According to Eq. (5.26), we 
can minimize the speckle effect by increasing either the laser beam 
divergence angle ( 0^) or the area of the receiver aperture (A R ), or both. 

However, the expression in Eq. (5.27) contains terms that are proportional 

2 ... 

to tan 9^. Consequently, the timing error will increase by increasing 


RMS TIMING ERROR IN PSEC 



NUMBER OF PHOTOELECTRONS (<N>) 


Fig. 5.5. RMS differential timing error versus received signal strength 
for pulse reflections from an infinite flat diffuse target 
using correlation estimator (z = 500 km, 9-p * 0.1 mrad , 

Of = 16.7 psec, Oj-, = 100 psec , * 0.5°, Xj = 0.355 urn, 

^2 = 1 .064 pm) . 
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0 T even though the 1 s increase. The only effective way to reduce the 
speckle-noise contribution is to increase the area of the receiver aper- 
ture. In fact, the increase in aperture size will also reduce the shot-noise 

1 

contribution because <N> is also proportional to A^. 

It is interesting to compare the performance of the correlation estima- 
tor with that of the optimal ML estimator. By using the fact that both the 
shot-noise and speckle induced signal fluctuations at the two frequencies 
are uncorrelated, the ML estimator of the two-color differential arrival 
time can be expressed as [Tsai and Gardner, 19851 


ML 


• f ?( t )dt 

J b 

-1 

+ 

' • f^(t)dt 

J s 

— <N 1 >“ 1 f g (t) + K^fgpU) 


-« <N 2 > _1 f s (t) + K2 1 f sp (t) 


-1 


(5.31) 


o? 


1 + 1 

<N X > g( <N L > , Kj) <N 2 > g( <N 2 > , K.j) 


where 


g«N.>,K 1 ) 


~ fdy 

/2rr -® 


exp ( — y 2 / 2 ) 


i-1,2 . (5.32) 


The timing error for the ML estimator consists only of first-order 
errors and so its dependence on the receiver bandwidth is only through o<, . 
Both Eqs . (5.29) and (5.32) are plotted versus the expected signal strength 
<Nj> = <N 2 > = <N> in Fig. 5.6 for comparison. 


II 
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NUMBER OF PHOTOELECTRONS (<N>) 


Fig. 5.6. Comparison of the JJMS differential timing errors using the corre- 
lation estimator ( T C0R^ anc * Maximum-Likelihood (ML) estimator 
( Tj^l) for pulse reflections from an infinite flat diffuse target, 
(z = 500 km, Of = 0.1 mrad, Of = 16.7 peec, $ = 0.5°, 

Xj = 0.355 urn, X 2 = 1*064 pm, Ar = 0.1 m^ , Kj ■ 9K£ = 10N 
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In the shot-noise-limited region (<N> « K^), the RMS timing errors of 

- 1/2 

the ML estimator have an <N> dependence as expected and are quite 

insensitive to the change in B . Whereas, the rms timing errors of the 

correlation estimator when B R >> B^ ( o^ ■ 50 psec, 100 psec) have an <N> * 

dependence due to the dominance of the second-order shot-noise error 

contribution. When B R is reduced ( o^ = 500 psec), the second-order 

contribution is less dominant and the overall timing error decreases as 
- 1/2 

<N> . In fact, the timing accuracy of the correlation estimator is 

comparable to that of the ML estimator when B R is chosen properly. 

In the speckle-limited region (<N> >> K^), however, we cannot gain much 
improvement by reducing the receiver bandwidth because speckle-noise 
contribution is much less sensitive to the change of B R . In this region, 
the ML estimator will perform consistently better than the correlation esti- 
mator . 

5 .4 Reflections from the Flat Cube-Corner Reflector (CCR) Arrays 

As the second example, we now consider a flat CCR array. The geometry 
and the orientation of the array is shown graphically in Fig. 5.7. The 
array has p CCR's in each row and q CCR's in each column. We shall assume 
the spacing d to be uniform between adjacent CCR'? in each row and in each 
column on the array, and the laser beam to be incident onto the target at an 
angle $ with respect to the target normal. The orientation of the array is 
such that the CCR's in each column are at the same range (height). This 
target configuration is similar to the CCR arrays on the Geodetic Earth 
Orbiting Satellite 1 and 2 (Geos-1 and -2) and the lunar arrays. 

For typical target ranges, the surface area of the CCR array is usually 
much smaller than the size of the laser footprint, and hence, we shall 
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(a) bottom view 


d sin $ 



Fig. 5.7. 


Geonetry of the flat CCR array 
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2 

assume the laser cross-section (|aj ) within the illuminated array to be 

uniform. Furthermore, we shall also assume the power reflectivities (8 ' s) 

m 

and the effective reflecting surface areas (A 's) of all the CCR's on the 

ID 

array to be equal. 

The total reflected electric field from a CCR array consists of the 
coherent additions of the field components due to the reflections from all 
the illuminated CCR's within the laser footprint. Due to the facts that the 
array is rough on the optical scale and that the range is continually 
changing, there will be slight variations in the apparent positions of the 
CCR's on the array. If the CCR array is m&croscopical ly flat and normal to 
the optical axis ($ = 0*), the distances traveled by the reflected field 
components back to the receiver plane may differ by several laser wave- 
lengths. As a consequence, the phase angles of the field components are 
random and mutually independent. This effect generates a partially deve- 
loped speckle field. The dephased components will introduce coherent inter- 
ference in the received signal power. On the other extreme, if the flat 
array is illuminated at a non-normal incidence angle so that the range 
difference between any two CCR's on the array is larger than the width of 
the transmitted laser pulse, then the received signal will consist of a 
train of non-overlapping pulses. In this case, the effect of coherent 
interference is negligible. In practice, the amount of interference will be 
somewhere between these two extremes. As expected, the partially developed 
speckle will introduce an additional error in the range measurements. 

The effect of coherency is given in Chapter 3 by the coherent 
interference factor y( p ,p ). The parameter which is most sensitive to this 
effect is the speckle SNR (K^). Because for satellite altitude ranging 
(Eq. (3.21)) is independent of the optical wavelengths, we shall let 
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Kj ■ Kj ■ K . By substituting Eq. (3.26) into Eq. (3.21) and using the above 
assumptions on the target parameters, K can be expressed as 
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m= 1 n 
m*n 




M 


M M 

l l 
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(5.33) 


where M“p*q is the total number of OCR's in the array, and ~ ijKj^) i- 3 

given by Eq. (3.29), and ji^ is the transverse coordinate vector of the mth OCR 
measured with respect to the center of the laser cross-section. Since the 
range delay due to beam wave curvature effects is small, Eq. (5.33) can be 
simplified to give 
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(5.34) 


For the target configuration given by Fig. 5.7, K can be expressed in 
terms of the OCR's spacing and incidence angle: 


(pq)‘ 


2p 2 1 (,-»)r ( " 6)2 . p 2 , - 


pq 


m“l 


(5.35) 


where 
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8 M d sin$ / (cOj) . (5.36) 

8 is the ratio between the range difference among OCR's on adjacent columns 
and the transmitted pulse width. When 8 * 0 (normal incidence), Eq . (5.35) 
reduces to 

K ■ pq / (pq-1) * M / (M-l) for 8 * 0. (5.37) 

In this case, K is always between 1 and 2 for M > 1. This low speckle SNR 
makes 3 peckle averaging impossible and is a consequence of the relatively 
small area of the OCR array. In other words, the speckle correlation area 
of the reflected field at the receiver plane is so large that only a few 
speckle cells fall within the receiver aperture. On the other hand, when 
d sin $ is much longer than the transmitted puise width (co^), Eq. (5.35) 
can be reduced to 

K ■ pq / ( p— 1 ) M M / (p-1), for 6 >> 1 . (5.38) 

According to Eq. (5.38), K can be very large if the OCR array consists of 
many widely spaced columns and only a few rows. In fact, if p * 1, then K ♦ " 
because in this case the reflected pulses from each CCR do not overlap in 
time at the receiver so there is no speckle induced interference. In 
reality, however, we would expect p and q to be equal and on the order of 5 
to 10 Con p -^uent ly , K will be somewhere between 1 and 10. Eq. (5.35) is 
plotted in Fig. 5.8 versus the horizontal spacing (d) of the CCR's within 
tne array for different incidence angles and for p=4 and q“3. 
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Fig. 5.8. Speckle SNR 00 versus relacive spacing of Che CCR array (d/co g ). 

Note chat Che relacive range difference between CCR's ac adjacent 
columns is (d/co g ) sin ^ where cc g is che point target pulse length, 
(z - 500 km, 0 T - 0.1 mrad. Of * 33 psec, - 33 psec, co g - 1.4 cm 
P - 4, q 
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The mean detected signal bandwidth can be calculated b/ using Eqs. 
(3.4), (3.24), (5.12) and (5.19). By neglecting beam curvature effects, the 
signal bandwidth has the form 


PT 


(1 - 


a d ) 
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a d } 


(5.39) 
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(5.40) 


(1-a^) is the bandwidth degradation factor and it is always between 0 and 

). For normal incidence, the reflected signal components from different 

OCR's will overlap in time almost perfectly. Consequently, is essentially 

zero, and the bandwidth of the reflected signal from the CCR array is the 

same as that from a point target. When the range difference between CCR's 

is on the same order of magnitude as the width of the point target response 

(co ), the signal components overlap partially and the total signal pulse is 

8 

broadened. In this case, we expect degradation in the detected signal band- 
width, and the amount of degradation depends on how much the range dif- 
ference between CCR's differs from ca . When the range difference is longer 

8 

than ca , the detected signal consists of a train of non-overlapping pulses 
8 

reflected from different CCR's. Although the detected signal is broadened 
in time, the bandwidth is still the same as that for a single pulse. 


n 
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The mean detected signal bandwidth for the configuation given in Fig. 
5.7 can be expressed as 
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S * 2a 2 
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(5.41) 


where 


8^ * d sin $ / 0 0^ 


(5.42) 


is the ratio between the range difference among CCR's on adjacent columns 

and the width of the point target response. Eq. (5.41) is plotted in Fig. 

5.9 versus the spacing of the array for different incidence angles $ and for 

p=4 and q=3. As predicted, the maximum bandwidth is obtained when $=0 or 

when the range difference (d sin $) is large compared with co . When 

g 

d sin <j) is on the same order of magnitude as co , however, the bandwidth is 

g 

degraded by almost a factor of 2. 

The rms bandwidth of the speckle variance (B<. p ) and the cross-bandwidth 
(B^) can be calculated by similar procedures. However, they both involve 
multiple summations, and numerical computations are usually required. 
Numerical results for Bg p and are plotted in Fig. 5.10 and 5.11, respec- 
tively, for the same parameters used in Fig. 5.9. Both B,. p and B^ have 
characteristics similar to that of the signal bandwidth B<,. 

The timing performance when ranging to the CCR array can now be 
calculated by using the above results for the speckle SNR and the various 
bandwidth parameters. Because the speckle noise contamination on the 
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Fig. 5.9. Normalized RMS detected signal bandwidth (B S /B PT ) versus relative 

spacing of the CCR array (d/cOg). Note that the relative range dif 
ference between CCR’s at adjacent columns is (d/c<Jg) sin where co 
is the point target pulse length, (z - 500 km, ■ 0.1 mrad, 

Of - 33 psec, 0 ^ » 33 psec, co 2 - 1.4 cm, p - 4, q - 3) 
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Fig. 5.10. Normalized RMS bandwidth of the speckle variance (B S p/Bp<j) verse* 
relative spacing of che CCR array (d/cc g ). Note that the relative 
range difference between CCR's at adjacent columns is (d/c<jg) sin 
where C 0 g is the point target pulse length, (z ■ 500 km, 

0 T - 0.1 rnrad, Of - 33 psec, ■ 33 psec, C 0 g - 1.4 cm, p - 4, q - 
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Fig. 5.11. Normalized RMS cross-bandwidth (D X /B PT ) versus relative spacing of the 
CCR array (d/co g ). Note that the relative range difference between 
CCR' s at adjacent columns is (d/cog) sin $ where co g is the point 
target pulse length, (z ■ 500 kn, 9-p ■ 0.1 mrad. Of ■ 33 psec, 

0 ^ - 33 psec, co g - 1.4 cm, p - 4, q - 3) 
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Fig. 5.12. RMS differencial tiding error versus received signal strength for 
pulse reflections from a CCR array. Note that 0^ is the relative 
range difference between CCR's at adjacent columns, (z ■ 500 km, 


> 0.1 mrad. Of * 33 psec, ov, 
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33 psec, co z - 1.4 cm. 
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Fig. 5.12. RMS differencial timing error versus received signal strength for 
pulse reflections from a CCS. array. Mote that 0^ is the relative 
range difference between OCR's at adjacent columns. (z ■ 500 km, 
9^ - 0.1 mrad, Of - 33 psec, ov, “ 33 psec, c a z - 1.4 cm, p - 4, 
q ■ 3) 
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In order to improve the ranging accuracy, very narrow transmitted pulses and 
receivers of very high resolution must be used. One important consequence 
is that the detected signal strength has very little effect on ranging 
accuracy. Thus, the design of two-color satellite ranging systems is domi- 
nated by the transmitted pulselength, pulse repetition rate and receiver 
bandwidth and not the detected signal strength. 

The above results can give us some insights on the timing accuracy when 
ranging to Lageos. Lageos is a 60-cm diameter sphere that is covered by 
426 OCR's. The spacing between adjacent OCR's is on the order of 10* arc- 
length with respect to the center of the satellite [Fitzmaurice efc al., 
1977]. This corresponds to about 5.2 cm. The corresponding range dif- 
ference between adjacent OCR's is approximately 4.5 mm, which is typically 
much smaller than the transmitted pulsewidth. Consequently, K is small, and 
the ranging accuracy to Lageos is also limited by the width of the detected 
signal pulse. 

5 . 5 Pulse Reflections from the Ocean 

The reflection of laser pulses from the ocean is due mainly to scat- 
tering by randomly distributed specular points on the ocean surface. These 
reflected pulses are broadened due to the differences in wave height (£(jj)) 
at various points within the laser footprint. The wave height at any given 
point on the ocean is the result of many wave components that have been 
generated by the wind in different regions and have propagated to the point 
of observation. Therefore, for footprint area which is large compared to 
the longest waves on the ocean, the central limit theorem can be used to 
argue that the surface profile of the ocean at the point of observation is 
Gaussian distributed. Under this assumption, the expected received signal 
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can be calculated by taking the expectation of Eq . (3.1) over the probabi- 
lity distribution of the Gaussian surface profile 5. In Fig. 5.13, the mean 
received signal by using the Gaussian ocean model is plotted for five dif- 
ferent laser beam divergence. As the beam divergence increases, the 
trailing edge of the pulse is broadened by the wavefront curvature effects, 
and the return pulse becomes asymmetric. 

The pulse shape computed using the Gaussian assumption will accurately 
predict the actual pulse shapes whenever the laser footprint is large com- 
pared to the ocean wave periods. However, when the footprint is small, the 
local probability distribution of surface profile within the footprint can 
be considerably different from a Gaussian. Sea wave records often reveal a 
sinusoidal or trochoidal profile for the dominant long-period waves. In 
these areas, the surface profile is more suitably modeled as a large ampli- 
tude sinusoid or trochoid wave with superimposed small amplitude disturban- 
ces (capillary waves) of short correlation length. For the one-dimensional 
sinusoidal wave model, the ocean surface profile can be represented by 
[Kinsman, 1965] 

£(_p) = A cos -jp x + <t>p + ^(jj) (5.45) 

where A is the amplitude of the sinusoidal wave, A is its period, is the 
phase angle, and (_p) is the small-scale disturbances. In Fig. 5.14, an 
expected received signal using this model is computed numerically for the 
case when the laser is pointed off nadir by 1°. The results reveal multiple 
glints (high frequency components) due to the specular reflections from suc- 
cessive wave crests and wave troughs. Crest reflections are responsible for 
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Fig. 5.13. Mean received pulse shapes for reflection from an ocean whose 
wave height is Gaussian distributed. The laser altimeter is 
pointed at nadir, and the significant wave height (SWH) is 1 m. 
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Altitude = 400 Km 

Transmitted Pulse Width = I cm ( 33ps ) 
Beam Divergence = 100/irad 
Wave Period = 10 m 

Wave Height = I m 
(Trough-Crest ) 

Nadir Angle <£ = l° 



J 
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Fig. 5.14. Mean received pulse shape for reflection from an ocean 
that has sinusoidal profile. The nadir angle is 1°. 
The wave height is 1 m (trough to crest) and the wave 
period is 10 m. The structure is caused by the 
reflections from successive wave crests and troughs. 
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Che sharp peaks in Che firsC half of Che pulse, while Crough reflecCions are 
responsible for Che sCrucCure in Che second half. The sharpness of rhe 
glinCs is relaCed Co Che beam divergence and Che ocean parameCers. 

Because of Che dynamic and random naCure of Che ocean, Che surface pro- 
file is noC a perfecc Gaussian or sinusoid. The Gaussian model and Che 
sinusoidal wave model can be viewed as Che Cwo exCremes. The acCual 
received pulse shapes aie likely Co lie somewhere beCween Che broad, smooch 
pulse and Che highly sCrucCured shapes. However, Che received pulses are 
expecCed Co exhibic fine scale SCrucCure due Co scrong specular reflecCions 
wichin Che laser fcoCprinC. This sCrucCure, which increasss Che signal 
bandwidch, is an imporCanC facCor in achieving picosecond ciming accuracies. 

To illusCraCe Che effecCs of surface profile on bandwidch (or Che 
sharpness of Che glinCs), Che MS bandwidch is calculaCed in Appendix B for a 
Gaussian disCribuCed profile under Che condicions ChaC Che profile correla- 
cion lengCh (L) is much less or greaCer Chan cne fooCprinC radius. The 
resulCs are 



A 1 ' 


A, / (a_ + 4 


2 . 2 \ 
o^/c ) 


■. 2. z tan 

(5.46) 


L << z Can 9^, 


where 


: T, ^ , 

^1 Y 1 
erfc(y.) 


(5.47) 



105 


e °e 2 2 2 ~ ^ 

2- Z -t ! n9 - : (2 V L + Can V 


L > z tan 9„ 


(5.48) 
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L << z tan 9„ 


and is the RMS wave height of the ocean. When L z tan 9 T> if the 
received pulse shapes are dominated by beam curvature (tan 9^ is large) or 
surface effects (o^ is large), Yj will be small.. If the curvature and sur- 
face effects are negligible, Yj will be large. On the other hand, when 
L « z tan 9_, y, will be large when o r is large. In these two limiting 
cases, Yj simplifies to 
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(5.49) 


2 

The other parameters such as the MS speckle bandwidth and the cross 

*\ 

bandwidth B“ can also be calculated using a similar procedure. The results 

A 

are tabulated in Table 5.1. 

To evaluate the expression for the timing ernr, we shall use the 
assumption that co^ << < z tan 9^. For large footprints the solution for 


TABLE 5.1 
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EXPRESSIONS OF THE QUANTIFIES WHICH ARE NEEDED FOR EVALUATING THE 
TIMING VARIANCE OF THE OCEAN-REFLECTED PULSES BY A CORRELATION RECEIVER 
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L << z tan 9^, is probably the most applicable. The MS timing error in this 
case is 
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(5.50) 


Since the second-order errors are not included, so Ineq. (5.50) only 
serves as the lower bound of the timing error. Because the ocean height is 
typically of the crder of tens of centimeters or greater, 0^/c will usually 
be a few nanosecond or more. Consequently, picosecond timing accuracy is 
probably not feasible when the footprint is large. 

When the laser footprint is small, the solution for L z tan 8^ is 
probably the most applicable, and the MS timing error in this case is 
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where 
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(5.52) 
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B q is the MS bandwidth of the signals that are reflected from the 
ocean. From Eqs. (5.51) and (5.52), we can see that when the footprint is 
small and the surface correlation length is Mrge, picosecond timing 
accuracy is possible if the bandwidth of the receiver is chosen properly so 
that it is not too large compared to the signal bandwidth. The effect of 
the receiver bandwidth on the timing accuracy for the ocean reflected pulses 
is similar to that in the diffuse target reflection situation. 

Consequently, the analysis on the effects of receiver bandwidth given in 
Section 5.3 can apply in the ocean reflection case. 

When the ocean surface correlation length is very large so that 
L » z tan 0 T , the ocean surface within the laser footprint is essentially 
flat. In this case, the MS timing error is 
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Notice that Eq. (5.53) is primarily the MS timing error for reflections 
from many specular points over a flat surface. This is expected because the 
ocean with very long surface correlation length behaves like a mirror within 
the laser footprint. 

In order to illustrate the effect of ocean surface correlation length 
(L) on the performance of the correlation estimator, the RMS timing error 
is plotted versus the normalized surface correlation length in Fig. 5.15 
for three different values of RMS ocean wave height ( o^) . The plotted 
results show that picosecond timing accuracy is achievable over regions 
where L is at least an order of magnitude larger than the laser footprint. 
When L is very large, the timing errors saturate at the lowest possible 
level . 

In Fig. 5.16, the RMS timing error is plotted versus the expected 
signal strength for three different sets of speckle SNR (K^'s) and when L is 
twenty times the footprint radius. The results show that the timing 
accuracy improves as the signal strength and speckle SNR increases. 

The timing accuracy of the ocean reflection using the sinusoidal wave 
model cannot be evaluated to closed form. Numerical results are plotted 
versus the expected signal strength in Fig. 5.17 for several values of 
speckle SNR = Kj - K. The system and ocean parameters correspond to the 


expected received pulse shape plotted in Fig. 5.14. Again, the timing 
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Fig. 5.16. RMS differential timing error versus received signal 
strength for the pulse reflections from a Gaussian 
ocean, (z = 400 km, 8-j« * 0.1 mrad, Of = 33 psec, 

0 ^ = 100 psec, X^ = 0.355 pm, Xo = 1.064 pm, = 0 , 

or = 3m). 
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PERFORMANCE OF CORRELATION RECEIVER 
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Fig. 5.17. RMS differential timing error versus received signal 
strength for the pulse reflections from a sinusoidal 
ocean. The mean received pulse shape is plotted in 
Fig. 5.14. 
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accuracy improves as the signal strength and speckle SNR increase. The 
system and ocean parameters in both Figs. 5.16 and 5.17 are deliberately 
chosen to be very close to each other in order to compare the timing accura- 
cies by using different ocean wave models. Although one was obtained 
through bandwidth analysis while the other through numerical computation, 
the results are very similar. According to these results, it is evident 
that the timing accuracy depends only on the bandwidth (or the sharp glint) 
of the reflected signal regardless of which ocean model that we assume. In 
both figures, picosecond timing accuracies appear to be feasible when the 
expected photocounts for both pulses approach a few thousand. The timing 
accuracy can be improved further by averaging the measurements from many 
shots. It should be emphasized that these results are only estimates of the 
expected timing accuracy. The actual receiver performance will depend on 
the local surface profile within the laser footprint. 

5 . 6 Discussion 

Pulsed two-color laser ranging systems can be used to determine the 
atmospheric delay by measuring the difference in propagation times between 
two optical pulses that are transmitted simultaneously at two different 
frequencies. However, proper design of the receiver timing algorithm is 
necessary in order to obtain picosecond differential timing accuracies which 
are required for most applications. In this chapter, the cross-correlation 
technique is used to estimate the differential propagation time between 
pulses which are reflected from CCR arrays and diffuse targets. The perfor- 
mance of the correlation technique is analyzed by considering both speckle 
and shot noise. The accuracy of this estimator is, in general, a function of 
the detected signal strength, speckle SNR, receiver bandwidth and the 
detected signal bandwidth. 



114 


Pulse reflections from various targets were analyzed in order to eva- 
luate quantitatively the performance of this estimator. For the flat diffuse 
targets, the timing accuracy is highly dependent on the receiver bandwidth 
and the characteristics of the time-resolved speckle. At low signal levels, 
it was shown that the timing performance of the correlation estimator is 
comparable to that of the optimal Maximum-likelihood (ML) estimator when the 
receiver bandwidth is chosen to match the bandwidth of the target reflected 
pulses. At high signal levels, however, time-resolved speckle places a fun- 
damental limit on the performance of the correlation estimator. For the 
cube corner reflector arrays, the timing performance is dominated by par- 
tially developed speckle so that the differential propagation time cannot be 
resolved to better than the pulaewidth of the received signals. For the 
ocean, the timing accuracy is highly dependent on the ocean surface correla- 
tion length and the bandwidth of the receiver. It was shown that picosecond 
timing accuracy is only feasible over regions where the surface correlation 
length is larger than the laser footprint radius. Based upon these analy- 
tical results, it appears the overall performance of two-color satellite 
laser ranging systems is dominated by the receiver bandwidth, transmitted 
laser pulsewidth, and pulse repetition rate, and not by the detected signal 


strength . 
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6. SYSTEM ALIGNMENT EFFECTS 

6 . 1 Introduction 

In the preceding chapter we described the performance of the correlation 
estimator under ideal situations. That is, we assumed that the two 
laser beams were perfectly aligned, that the beam divergence angles asso- 
ciated with the two operating frequencies were identical, and that there was 
no displacement between the transmitter and the receiver. In actual prac- 
tice, however, not all the above assumptions will be satisfied, and as the 
result the timing accuracy we obtain may be degraded. In this chapter we 
investigate the system alignment effects on the performance of a two-color 
laser ranging system. 

6.2 Correlation Coefficient and Mean Square Error 

One of the basic assumptions used in applying the correlation timing 
estimator is that the return signals at two frequencies have similar pulse 
shapes. In practice, however, the return pulses may not always be iden- 
tical. This could be due to misalignment at the transmitter or due to 
atmospheric turbulence. Any decorrelation of the pulse shape will, in 
general, degrade the performance of the correlation estimator. 

The correlation coefficient between the received signals at the two 
frequencies can be defined as 


p 12 <T) 
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( 6 . 1 ) 
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When the two signals are identical in shape, has a maximum value of one 
at T - x , where x q is the difference between the propagation times of the 
two signals, p^ 2 is in general random due to the random nature of the 
received signals. The expected value of p^ is approximately given by 


<P 12 (t)> 2 


/ < S^t) S 2 (t+T) > dt 
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[/ dtj / clt 2 < S 2 (tj) S 2 (t 2 +t) >] 
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( 6 . 2 ) 


Using the fact that the shot noise and speckle are uncorrelated at the two 
frequencies, the expected value of p^ 2 can be simplified to 
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When the signals are not identical, Eqs . (5.8) and (5.9) become 
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where fg and fgp are Che normalized received signal and speckle variance 
1 i 

at wavelength A^ , respectively. By substituting Eqs. (6.4) and (6.5) into 
Eq. (6.3) and using Parseval's Theorem, we can express <Pj 9 (t)> at t“T q in 
terms of the frequency spectrum of the mean received signals as 
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When the mean received signals at the two frequencies are different, 
the MS timing error of the correlation receiver becomes 



( 6 . 8 ) 


The numerator of Eq. (6.8) can be expressed in terms of the Fourier 
transforms of the expected signals and the speckle variances as 
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where 



/ dt f sp (t) e~ JUit 


i-1,2 


( 6 . 10 ) 


The first term in Eq. (6.10) equals the square of the bias. The next 
four terms are the first-order terms involving the cross correlations of the 
mean of and the fluctuating portion of S^, and vice versa. The last four 

terms are the second-order terms involving the cross correlations of the 
fluctuating portions of the two signals. 

6.3 Geometry 

The geometry of the ranging system is illustrated in Fig. 6.1. The 
laser beams at two different wavelengths are transmitted from the system 
simultaneously. Laser beam 1 is incident on the target surface at an angle 
$, and beam 2 is incident on the target surface at an angle $ + A$. A$ is 
the difference in incidence angles between the two beams due to misalign- 
ment. Throughout the analysis, the incident angles are assumed to be small 
so that the shadowing effects can be neglected. The expressions for the 
means and covariance of the received signals involve 2-D integrations over 
the planes that are transverse to the direction of propagations. Therefore, 
we need to determine the apparent distances (z.) and the apparent surface 
profiles (5(jO) in terms of z and £(_p) . From simple geometrical con- 
siderations, we have 
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* Z / cos ( $ + (i-l)A$) , i*l ,2, (6.11) 

x. = x cos( + (i-l)A$) - £(j>) sin( $ + (i-l)A$) (6.12) 

=« x cos($ + (i-l)A$) , i* 1,2 

y£ * y . i*l,2, (6.13) 

5(j>£) * x sin( $ + (i-l)A^) ♦ 5(j)) cos($ ♦ (i-l)A$), i*l,2. (6.14) 


The laser intensity footprints at the two wavelengths are assumed to be 
Gaussian, 
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where 


* z tan(d T + (i-l)A8 T ) / cos($ + (i-l)A$) . (6.16) 

- 1/2 

o^ is the footprint radijs measured at the e point on the _p^ plane. In 
Eq. (6.15), we have assumed that laser beam 1 has a divergence angle of 
9^ and beam 2 has a divergence angle ( 8^ + A8^.) . A9^, is the difference in 

divergence angles between the two beams. 

The power spectra of the normalized received signals can be calculated 
by using Eq. (3.24), Eq . (6.7), Eqs. (6.11) through (6.15) as 
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where = (x a> y a ) and ^ - (x b ,y b ). 

For the Gaussian target surface profile with RMS surface roughness a 
the ensemble average in (6.17) becomes 
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where R^( p^,^) is the normalized autocorrelation function of the surface 
profile. If the correlation length of the target surface (L) is large 
compared to the laser footprint radius, R^J^.J^) can be approximated as 
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Substituting Eq. (6.19) into Eq . (6.11) and carrying out the integration, we 
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have 
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The cross spectrum between the two signals (6 $ ) can be calculated 

S 1 S 2 

using a similar approach with slight modifications of the expected phase 
function. Due to the difference in incidence angles at the two propagation 
paths, the laser foortprints do not overlap completely, and R^( _P jt _P^) must 
be modified so that 
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After carrying out the integration, we have 
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where 
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The expected correlation coefficient can now be evaluated by substi- 
tuting Eqs. (6.20) and (6.23) into Eq. (6.6). Numerical results are plotted 
versus different system parameters in Figs. 6.2 and 6.3. In these results, we 
let <p - 0 because it is the most common viewing configuration in satellite 
laser ranging. Similar results vrould be obtained when p is nonzero but small. 
In Fig. 6.2, <p^> is plotted against 48^/8^ for different values of misalign- 
ment angles 4$ and for a target surface correlation length which is on the 
same order of magnitude as the smaller footprint. We can see that the maximum 
correlation can be obtained when both 48^ and 4$ are zero; that is, when the 
two beams have the same divergence and they are aligned so that the laser 
footprint intensities are identical. This maximum value, however, is less 
than unity because of the decorrelation of the signals due to shot noise and 
spccxle (refer to Eq. (6.6)). When the two beams are misaligned, the 
footprints overlap only partially, which results in differences in received 
pulse shapes and a decorrelation of the two c'gnals. The peak correlation at 
each curve corresponds to the situation in which the largest possible overlap 
area is attained. 

In Fig. 6.3 * s Pl° tte d versus the target surface correlation 

length and for 48^ =* 0 (equal divergence of the two beams). When the magni- 
tude of the surface correlation length is on the same order as the size of the 
footprint, the degree of correlation is roughly proportional to the overlap 
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Fig. 6.2. Mean correlation coefficient versus difference in divergence 
angles . 
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SURFACE CORRELATION LENGTH IN METER 


Fig. 6.3. Mean correlation coefficient versus target surface 
correlation length. 
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area. Consequently, the top curve, which corresponds to the largest overlap 
area, has the highest correlation. When the surface correlation length is 
very large (flat surface), the intensities reflected from different regions of 
the surface are identical. In this case, high correlation can always be 
achieved regardless of how well the two footprints overlap. 

In Chapter 5, it was shown that the timing performance of the correla- 
tion algorithm depends on the ratio of the receiver bandwidth (B^) to the 
expected signal bandwidth (Bg). Since this algorithm depends on how well 
the two signals are correlated, we could expect the correlation coefficient 
to exhibit the same dependence. We shall illustrate this point by means of 
a simple example. Let us assume that the received signals have Gaussian 
pulse shapes and that there are no system alignment errors. In this case, 
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where 
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Og is the RMS pulse width of the received signal at both wavelengths and 
0^ is pulse broadening due to the target range spread. It follows from Eq. 
(6.25) that Bg and are 
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and 
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respectively. By substituting Eq. (6.25) into (6.16), we have 
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When ranging to diffuse targets at satellite altitudes, is much 
larger than <N^>. If we let <Nj> = <^> = <N> , Eq . (6.26) can be simplified 
to 
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Notice that the receiver bandwidth B R is always greater than the 
detected signal bandwidth Bg. According to Eq. (6.27), B R must not be too 
large compared to Bg in order to obtain high correlation between the two 

t 

signals. This is the same conclusion reached in Chapter 5. 

The numerical results for the MS timing error are plotted versus the 
signal strength in Figs. 6.4 and 6.5 for different ranging parameters. 

For simplicity, we have assumed the signal strengths at the two wavelengths 
are equal. In Fig. 6.4 we assume the two beams are aligned but their 
divergence angles are different. The bottom curve (49^, = 0) corresponds to 
the case of complete overlap of the laser footprints (no system error). As 



oruou z-« TDOxj^m nz-z«H win 


5 

3 


Z = 400 km 
cr f = 16.7 psec 
<r h = 100 psec 
# T = 0.1 mrad 
L = 50 m 
a£ = 3m 
X, = 0.355/im 
X 2 = 1.064 ^im 
A 4 >= 0 
K, = 9K 2 = I0 5 
£ = 0 


102 i03 104 i 0 5 i 0 6 10 7 i00 i0 q 

number of photoelectrons 


RMS differential timing error versus expected signal strength when 
the divergence angles of the two transmitted laser beams are dif 
ferent . The ranging system is pointed at nadir. 


expected, the timing error increases wi b increasing This additional 

error is partly due -> the slight difference in signal intensities introduced 
by different footprint areas, and partly due to the broadening of one of the 
received signals. However, this additioual error is usually small. 

In Fig. 6.5 we assume the two beams have equal divergences but they 
are misaligned. The remaining parameters are the same as those used in 
Fig. 6.4. As we can see, the additioual error due to misalignment is con- 
siderably larger than that due to differences in divergence. The middle 
curve corresponds to the case of partial overlap of the beams, and the addi- 
tional error due to this slight misalignment is not very severe. The top 
curve corresponds to the case at which the two beams illuminate different 
parts of the target surface. Since the surface correlation length is 
finite, the reflected signal intensities at the two wavelengths are con- 
siderably different. Consequently, the performance of the correlation 
receiver in this case is degraded considerably. 

6 , 4 Return Signals for Coaxial and Noncoaxial Systems with Central 

Obstruction 

For the laser ranging systems, the recieved signal is given by the link 
equation 


<»> ' h 0 r « T a \ * 


(6.28) 


where Q is the average transmitted pulse energy, n is the quantum efficiency 
of the receiving system, hf is the signal photon energy, T a is the two-way 
atmospheric transmittance, 0° is the average power reflectivity of the 
target, A^ is the receiver aperture area, and z is the range distance. In 
most cases, however, not all the light collected by the telescope is 
detected by the detector. Signal power is normally lost by two mechanisms. 
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First, the overlap of the transmitted beam with the receiver's f ield-of-view 
v.FOV) is often Incomplete as in noncoaxial systems or in systems with sha- 
dowing of the main mirror by a central obstruction, e.g., a secondary 
mirror. Second, because the telescope is focused at infinity, light 
backscattered from targets at email and medium distances is not focused in 

the focal plane. Both effects give rise to a deviation of the detector 
-2 

power from the z dependence as given by the link equation. 

To study the return signals for coaxial and noncoaxial systems, we 
assume a target with uniform reflectivity and negligible range spread. The 
analysis here closely follows that of [Harms et al., 1978, 1979]. 

The geometry of the receiving optics is given by Fig. 6.6. Let 
r^ and r^ be the radii of the primary mirror and central obstruction, 
respectively. The detector is situated on the optical axis, at a distance D 
behind the primary lens, while the obstruction or secondary mirror of the 
telescope is located at a distance in front of the lens. The primary 

lens forms an image of the target on the image plane. For the target at 
distance z, the image is located at distance I(z) behind the lens, where 
I(z) is related to z by the thin lens formula 

1/I(z) + 1/z - 1/f , (6.29) 

and f is the focal length of the primary lens. 

The light backscattered from a point on the target is focused onto 
the image plane resulting in an irradiance dS(£^,z) at point _p^ , where 
£^= 2 ^ I(z)/z. A pencil of light backscattered to a surface element dA 
around £^ will illuminate an annulus in the detector plane around _p^ with an 


irradiance of 



Fig. 6.6. Geometry of the receiving optics. 
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(6.30) 


where S(_p^,z) and S(jk,z) are the irradiances in the detector and image 
plane, respectively. The size of the annulus can be calculated by using 
geometric optics. The ray tracing of the reflected light is given by 
Fig. 6.7. From this figure, ^ is related to jk by 


-2d = im-2i 


(6.31) 


and the outer radius of the annulus r is related to r by 
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The total irradiance at the detector plane is provided by all the pencils of 
rays for which the principal rays intersect the detector plane within a 
circle around with radius r £ excluding those rays within a circle with 
radius r c r^ / r^. The corresponding marginal points in the image plane 
form an annulus with an outer radius r^ given by 


r - r iii! 
m c D 


(6.33) 


and an inner radius equal to r r./r . 
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Fig. 6.7. 


Shadowinu in the detector plane due to central obstruction. 
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In the noncoaxial systems, the position of the transmitted beam rela- 
tive to the receiver axis varies with distance z. Fig. 6.8 shows the 
geometrical arrangements of the transmitter and receiver in a noncoaxial 
ranging system. The detector defines the optical (z) axis of the receiver, 
the y axis points toward the transmitter, and the x axis is perpendicular to 
both. The transmitted beam and receiver axes are inclined toward eacn other 
by a small angle 6 with a separation d at z = 0. Therefore, * - he 
y-coordinate of the beam center at any distance z is given by 

a(z)*d-6z (6.34) 

o 

Because of the displacement, the total irradiance distributions on both 
the image plane and the detector plane are no longer symmetric around the 
receiver axis. Instead, they are symmetric around a parallel axis whose 
y-coordinate is displaced from the receiver axis by the amount 

a( z) = (d - 6z) . (6.35) 

z 


By assuming a Gaussian footprint with radius z tan 0^ at the 

- 1/2 

e point, and using Eq. (6.35), the irradiance at the image plane is 
given by 
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where 


z tan 0„ 


I(z) 


15 I(z)tan 9 


(6.37) 


— i = ^ x i’ ^ » and P 0 c ^ e average transmitted power. 

By using Eqs. (6.30) and (6.36), we can express the irradiance at the 
detector plane as 
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where A(r , P.)/A(r — , — P.) denotes the integration region - an annulus 
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By making changes of variable, we can write 
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The average detected power is given by 

2ir r, 

P(z) = / d_p^ S(_p^ ,z) - J d$ / p d dp d SU^.z) (6.40) 

where r^ is the radius of the detector. By using the assumption that 
d << z r^ / D, Eq. (6.40) can be evaluated to give 
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where I is the modified Bessel function of the first kind and zero order, 
o 

Again, by making changes of variables, we have 
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The expected signal strength is then 
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If the radius of the detector r d is very large, i.e., r d » z tan © T , 
Eq. (6.43) can be simplified to 
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(6.44) 
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In Eq. (6.44), the term (r - r.) includes the effect of central obstruc- 

p D 

tion; the exponential term includes the effect of separation between 
transmitter and receiver. It can be seen that the familiar link equation 
can be deduced for the coaxial system when the detector size is large, 
except that in this case the effective receiver aperture is being reduced by 
the obstruction. For the noncoaxial systems, on the other hand, the 
received signal strength can be considerably different from that predicted 
by the link equation. 

6.5 Numerical Results 

For the coaxial systems, the change in the return signal strength is 

mainly due to central obstruction and finite detector size. In Fig. 6.9 

the received signal strength is plotted versus altitude for the coaxial 

system without central obstruction and for different detector radii. The 

telescope is assumed to be focused at infinity. When the detector is very 

large, the result reduces to the familiar link equation so that there is a 

z -2 dependence. For a finite detector size, signals at close distances are 

partially detected because the target is not focused on the detector. In 

Fig. 6.10, we show the results for the system with central obstruction. 

For large detector size, the return signal is uniformly reduced in propor- 

-2 

tion to the smaller effective receiver area. However, the z dependence is 
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Fig. 6.10. Expected signal photocount versus range for 
central obstruction. Telescope is focussed 


the coaxial system with 
at infinity. 
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maintained because all the energy collected by the telescope Is detected 
(Eq. (6.44) with d * 5=0). For small detectors, on the other hand, the 
signals reflected from targets at close distances are reduced appreciably 
because the obstruction causes a shaded region in the detector plane. This 
shaded region is directly related to the size of the obstruction. 

In the noncoaxial systems, the change in return signal strength is 
mainly due to the relative orientation between the transmitter and receiver. 
Fig. 6.11 shows the received signal strength as a function of distance for 
a system separation (d) of 9 inches and for different inclination angles 5. 
In this figure, we do not include the effect of obscuration. For com- 
parison, the results from the standard link equation and the signal levels 
of a corresponding coaxial system are also given. It can be seen from the 
figure that even for the modest separation of the transmitter and receiver, 
only the signals of limited range region are detected by the receiver. The 
strongest signals are detected from the distances where the laser footprints 
and the receiver FOV's completely overlap each other. The signal levels 
fall off rapidly as the overlapping areas decrease. The distance at which 
a peak signal can be detected decreases with increasing inclination angle 5. 

In Fig. 6.12 the signal strength is plotted versus distance for the 
noncoaxial system with central obstructon. The parameter used in this 
figure is the same as that in Fig. 6.11. We can see the central obstruc- 
tion reduces the signal level uniformly at large ranges. In addition, the 
return signals from short distances are reduced in strength due to the sha- 
dowing in the detector. However, the effect is relatively minor compared 
to the effect of system separation. 
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Fig. 6.11. Expected signal photocount versus range for the noncoaxial system 
without central obstruction. Telescope is focussed at infinity. 

The transmitter and receiver are separated by 9 inches. The remain- 
ing parameters correspond to those given in Fig. 6.9. 
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Fig. 6.12. Expected signal photocount versus range for the noncoaxial system 
with central obstruction. Telescope is focussed at infinity. The 
transmitter and receiver are separated by 9 inches. The remaining 
parameters correspond to those given in Figure 6.10. 
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In Fig. 6.13, the received signal strength is plotted versus the 
displacement d for the system with central obstruction and with parallel 
axes (6 = 0). This system, although ideally matched at infinity (complete 
overlap), suffers from a strong reduction of the near range signals. This 
reduction of the near range signals is not due to the obstruction, but is 
the result of incomplete overlap between the transmitter and receiver at 
close distances. Consequently, this system is suitable for satellite 
distance ranging, 'nit not so for airborne ranging. 

6.6 Conclusion 

T n this chapter the system alignment effects on the performance of 
two-color laser ranging systems were investigated. When the transmitted 
laser be -in at two wavelengths have different beam divergences and/or are 
misaligned, the correlation between the two signals is reduced considerably. 
Consequently, the timing performance of the correlation receiver is 
degraded. The degradation is particularly severe when the two beams are 
misaligned. It was also shown that the expected correlation coefficient 
depends on the expected signal strength and the relative magnitude between 
the receiver bandwidth and the bandwidth of the mean signal. The best 
correlation can be obtained when there is no system error, when the signal 
strength is high, and when the receiver bandwidth is only slightly larger 
than the signal bandwidth. 

The effects of the central obstruction on return signal strengths were 
also investigated. In coaxial systems, the central obstruction is more 
important to signals from short distances because it causes a shaded region 
in the detector plane and results in appreciable reduction in signal level. 
For large distances the obstruction gives rise to a uniform reduction of the 
signal level. In noncoaxial systems, the suppression of return signals is 
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Fig. 6.13. Expected signal photocount versus range for the noncoaxial ranging 
system with parallel transmitter and receiver axes (6 = 0). 
Telescope is focussed at infinity. The remaining parameters 
correspond to those given in Figure 6.10. 
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primarily due to the incomplete overlap of the transmitted beam with the 
receiver’s f ield-of-view. Noncoaxial systems with parallel axes suffer from 
considerable loss in signal return at small to medium distances and approach 
the coaxial signal for long distance. The obstruction effect of the 
receiver appears to be less significant as compared to the effects of the 
geometrical arrangements in the noncoaxial systems. The advantage of using 
the noncoaxial systems is that the unwanted signals (such as background 
noise) can be compressed by choosing the system separation and inclination 
angle carefully so that complete overlap between the laser footprint and 
the receiver FOV occurs only at the target plane. 
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7. HORIZONTAL, PATH RANGING EXPERIMENTS 


7 . 1 Introduction 

Parallel with the theoretical research on the pulsed two-color laser 
ranging, an experimental program has been administrated by NASA-Goddard 
Space Flight Center (GSFC) to determine the performance of the available 
system hardware and identify possible hardware limitations in the past five 
years. By using a mode-locked ND:YAG laser transmitter and a 
photomultiplier-based receiving system, Abshire [1982] first demonstrated 
that differential timing accuracies of a few picosecond was achievable when 
ranging to a single cube corner reflector over a horizontal path of 500 m. 

In March of 1983, we conducted a single-color horizontal path experiment at 
GSFC to study the speckle effects and timing accuracies by ranging to a flat 
diffuse target and a ocean-simulated target [Im, et al, 1983]. Our results 
showed that single shot accuracies of 90 to 40 psec, respectively, were 
acheiveable for these two targets. In September of 1983, the first two- 
color airborne experiment was conducted by ranging to the ocean [Abshire, 
et al., 1984, Tsai and Gardner, 1984], The results indicated high correla- 
tion between the returns at the two wavelengths, and the single shot 
accuracy approaching 40 psec was achievable. 

The above experimental results indicate that the two-color laser 
ranging technique is feasible. Yet, the ranging accuracy at that time was 
limited by the relatively long impulse response time (500 psec) of the 
system electronics. For the technique to apply to satellite laser ranging, 
a receiver of much higher time resolution is desirable. In late 1983, a 
dual-channel streak camera of 5-psec time resolution was interfaced into the 
ranging system by the NASA personnels. In February and April of 1985, a 
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series of horizontal path experiments were conducted using the new system. 
The main objective of these experiments was to verify the theory behind the 
cross-correlation differential timing algorithm (Chapter 5). In this 
Chapter, the experimental measurements will be presented and compared with 
the theory. 

7.2 Review of Theory 

In a two-color pulsed laser ranging system, the laser pulses of two 
different wavelengths are transmitted simultaneously. Due to the difference 
in refractive indices of the atmosphere at different optical frequen- 
cies, the two pulses will arrive at the target surface at different times. 

If the laser intensity cross-sections at the two wavelengths are approxima- 
tely equal, the mean reflected pulse shapes detected by the receiver will be 
very similar. Therefore, the difference between the roundtrip propagation 
times of the two signals can be estimated by cross-correlating the detected 
pulses and searching for the peak of the correlation function. 

Mathematically, the cross-correlation estimator can be written as [Im and 
Gardner, 1985] 

= arg max I Rj 2 ^ ^ * (7.1) 

t 

Sj and are the detected signals at wavelengths X^ and » respectively. 

R, 2 (t) is the cross-correlation function of these two signals. T is the 

differential propagation time, and T is its estimate using the correla- 
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tion algorithm. Arg max (f(x)) is used to denote the argument x that maxi- 
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The performance of the correlation estimator was derived in Chapter 5. 
This estimator is unbiased, and the MS timing error is highly dependent on 
the bandwidth of the mean received signal, receiver bandwidth, expected 
signal strength and the time-resolved speckle. 

The bandwidth of the received signal is inversely related to the time 
extent of the received signal, which in turn is related to the ranging 
geometry and the reflecting properties of the target surface. When ranging 
to a diffuse target, the reflected signals are broadened mainly due to the 
range spread of the target. When ranging to a specular target such as a 
cube corner reflector (CCR) array, the reflected signals contain sharp 
glints which characterize the point target response of the ranging system. 
The signal statistics for these two target configurations were given in 
Chapter 3. For the experimental work reported here, we are primarily 
interested in the pulse reflections from a tilted diffuse flat target board 
with its surface area larger than the size of the laser footprint, ar.d from 
a tilted CCR array with its surface area much smaller than the size of the 
laser footprint. The latter configuration is of particular interest 
because CCR arrays have been deployed on satellites (LAGEOS and GEOS) and 
the moon and are routinely used as targets for geodetic measurements. By 
assuming Gaussian models for the receiver impulse response and the 
transmitted signal intensity, the expected received signals for these two 
target configurations can be expressed as 


153 


<S(t)> 


- t 2 llal 


<N> 


/57 


<N> 


M 


/Tit a m-1 

g 


exp 


for tilted diffuse flat target 


[t - 2 C(d)/c]‘ 
— <n 

2 o 2 


(7.2) 

for tilted CCR array 


where 


, 2 2 , 1/2 

<0 g 4 °r> 


,„2 . „2 2 . 1/2 

(<J f + °h + a r ) 


(7.3) 


is the RMS pulse width of the detected singal that has been reflected from a 

tilted diffuse flat target, and is the pulse broadening due to the range 

spread of the target. The remaining quantities were defined in Chapter 5. 

According to Eq. (7.2), the mean received signal reflected from a tilted 

diffuse flat target is Gaussian shaped, while the mean received signal 

reflected from a CCR array consists of a train of narrow pulses with widths 

equal to the point target response of the receiver. Each pulse corresponds 

to reflections from one CCR. When C = 0, the individual pulses merge to 

form a single pulse. When 0 < 5 < co , the pulses overlap partially which 

g 

results in broadening of the overall pulse shape and a reduction in signal 
bandwidth. When £ > co^, the pulses do not overlap so that the signal band- 
width is equivalent to that of a single pulse. The RMS pulse width of the 
reflected signal from the CCR array can be calculated from Eq. (7.2) as 
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where £ is the average range of the CCR array 
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2 

The power spectrum of the mean recieved signal (|$g(u)| ) is par- 
ticularly useful in determining the mean signal bandwidth. It can also pro- 
vide information on the photocount of the mean pulse measurements provided 

2 

the system is properly calibrated. |$g(w)| is defined as 


|'t> s (w) 


|2 _ 


/ S(t) 




dt 


(7.6) 


For the flat diffuse target and the CCR array, Eq . (7.5) can be evaluated 
explicitly using Eqs. (7.2) and (7.3) to give 
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where 


]G(id) | 2 - |H(u) | 2 |F(u) | 2 


(7.8) 


G(u>) is the Fourier transform of the system point target response, H(u>) is 
the transfer function of the receiver, and F(u>) is the Fourier transform of 
the transmitted signal. It follows from Eq. (7.7) that the MS bandwidth of 
the mean received signal (Bg) is 
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tilted CCR array 


where (l-cx^) is the bandwidth degradation factor derived in Eq. (5.40). It 

takes on values between 0 and 1. When the pulse reflections from different 

OCR's completely overlap (5 = 0) or when they are completely separated 

(5 > co ) in time, is zero and the bandwidth is equal to that of a single 

pulse. When the pulse train overlaps partially (0 < £ < co ), the mean 

S 

signal bandwidth is reduced because (1- a^) is less than 1. 

When the range spread of the target is larger than the correlation 
length of the speckle induced fluctuations, the time-resolved speckle will 
cause small-scale fluctuations within the detected pulse and distort its 
shape. In this case, it may be difficult to extract the mean signal from 
the measured data accurately. We can, however, obtain a biased estimate of 
the mean signal power spectrum by averaging the power spectra of the 
measured pulses. This technique is particularly useful in situations where 
the mean signal is not known a priori such as ranging to dynamic targets. 

The expected power spectrum of the received signal is defined as 


E [ | S( u>) | 2 ] = < | S ( u») | 2 > 


(7.10) 
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where S(u>) is the Fourier tranform of the S(t). After expansion, we have 


< |s(u>) | 2 > - 


• " - joKtj-t,) 

/ / dt 1 dt 2 <S(tj) S(t 2 )> e 


(7.11) 
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where C c is given by Eq. (3.23) and is the covariance of S(t). From Eqs. 

u 

(3.22) and (3.23), it follows that 


< | S ( o») 1 2 > = <N>|H(u))| 2 + <N> 2 K _1 |G(w) I 2 + |* s (m) | 2 


(7.12) 


where |h| 2 and | G | 2 are the power spectrum of h(t) and g(t), respectively, 

2 

and K is the speckle SNR. As can be seen from Eq. (7.12), <|s| > is a noisy 
o 

estimate of |$| . The first term is due to shot noise, and the second term 

is due to speckle. These two terms are most significant at high frequen- 

2 

cies. To reliably estimate | * s | , the signal must be large compared to the 
shot noise and speckle level. After combining Eqs. (7.7) and (7.12), we 
have 


<N> | H( m) | 2 + <N> 2 |G(u>) 


2 2 

, - u» a 

I + e 

K 


tilted diffuse flat target 


< | S( co) | 2 > = 


<N> | H( u>) | 2 + <N> 2 |G(u») | 2 


tilted CCR array 


1 1 r v 

7 + “7 L l cos 

* M m-1 n*l 


M M 


— (5(0 ) - 5(p )) 

c “in — n 


(7.13) 


158 


For the diffuse target, l*g|2 can be estimated accurately at the frequency 
range 


1 1/2 
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We can see from Ineq. (7.14) that <|S((»)| s can estimate |$ | accurately 

2 ) 

when <N> and K are large, and when the received signals have narrow pulse 

2 

widths. For the CCR array, accurate estimate of | <5g | can only be achieved 
when the array surface is normal to the optical axis such that 


5(p) - 
— (11 


5(0 


V m,n, 


(7.15) 


and when both <N> and K are large. When the array is tilted, there will 
always be some frequency components 


01 


kirc 

U(p > - 5<p. 
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(7.16) 


that fall below the shot noise level and cannot be estimated. 

The various parameters mentioned in this section will be compared with 
the experimental measurements in subsequent sections. 


159 


7 . 3 Experimental Configuration 

The major elements comprising the two-color laser ranging system used 
during the experiments were a mode-locked, frequency-triple ND:YAG laser 
transmitter and an image-intensified streak camera-based receiver. 

The laser transmitter assembly is given in Fig. 7.1. When the laser is 
triggered, its flashlamps are fired, and the 1.064-um radiation rapidly 
increases within the oscillator cavity. This cavity contains both a dye 
cell to force mode-locking, and an etalon to control the width of the 
oscillator pulse. The optical output from the oscillator contains 6 to 12 
of these pulses with a pulse separation time of 7 nsec. The Pockels cell is 
triggered by the oscillator output pulse train. It and the output polarizer 
pass only a single output pulse and attenuate the remaining input pulses by 
20 times. The energy of the selected pulse is then increased by ten times 
in the amplifier assembly. This output pulse is then input to the frequency 
doubler and the tripler subsequently which convert part of its energy to 
0.532 and 0.355 pm. The nominal laser output is 6 mJ at 1.064 pm, 3 mJ at 
0.532 pm, and 1 mJ at 355 pm. The three output pulses are coincident in 
time and exit the transmitter simultaneously. 

The streak-camera-based receiver uses the swept-image converter to pro- 
vide psec response time, which is about 15 times faster than the transmitted 
pulse widths. The operation of the linear scan streak camera is shown in 
Fig. 7.2. The incident optical pulse illuminates the p'lotocathode and frees 
phocoelectrons from its rear surface. These photoelectrons are accelerated 
rapidly by a mesh electrode and are deflected vertically by a fast electrical 
sweep. This results in a time-to-space mapping of the electron stream. The 
resulting electron distribution impinges on the microchannel plate, which 
preserves the spatial distribution and amplifies it. The amplified electron 
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bundle exits the rear of the plate and is accelerated into a phosphor 
screen. This produces a weak optical image whose spatial intensity distri- 
bution is proportional to the temporal intensity distribution of the illumi- 
nating pulse. An intensified video camera reads out this image, and a 
digitizer-microprocessor system converts it back to an intensity versus time 
profile . 

The ranging system configuration used during the CCR array test is 
shown in Fig. 7.3. During the diffuse target test a different transmitter- 
telescope orientation was used. This new subsystem is shown inside the 
circle, which replaces the subsystem within the dashed box. However, the 
principal operation of these two systems are identical. The system parame- 
ters are listed in Table 7.1. 

The systems were operated in the following manner. The 30 psec (FWHM) 
pulses at both 0.532 and 0.355 pm from the YAG laser first entered the 
transmit beam alignment system to eliminate the angular offset due to the 
nonlinear action of the angle-matched KD P freouency doubling and tripling 
crystals inside the laser. Four percent of the aiigned fignal energy was 
then reflected by a beam splitter into the reference path to start the 
timing system and to determine the system delay. The remaining signals at 
the two wavelengths were directed to the target along a horizontal path. 

The optical returns from the targets were collected by the receiver 
telescope. They were reimaged by an eyepiece, which relayed the focal point 
to the start of the optical delay system. At this point, a beamsplitter 
reflected part of the optical signals into a photomultiplier, whose output 
was used to trigger the streak camera and to stop the timing system. The 
remainder of the reflected energy at each wavelength entered the 40 nsec 
imaging delay line. This additional delay allowed enough time for the 
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Fig. 7.3. The pulsed, two-color laser ranging systems used during the 
CCR array and diffuse target experiments. 
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TABLE 7.1. 


SYSTEM PARAMETERS FOR THE TWO-COLOR LASER RANGING SYSTEMS 
USED DURING THE HORIZONTAL PATH EXPERIMENTS 


System Parameters 
Horizontal Path experiments 


Laser 


Quantel YG40 - Mode Locked ND:YAG 
nominal 30 psec pulse (FWHM) 

0.3 raJ at 0.532 pm, 0.5 mJ at 0.355 pm 


Telescopes 


CCR Array Test - 460 cm are^, 91 cm focal length 
Diffuse Target Test - 891 cm area, 391 cm focal 

length 


Mirrors 


Double Dielectric 

maximum reflections at both 0.355 and 0.532 pm 


Photomultiplier 
Optical Delay 
Streak Camera 


TIU 


Computer 


Hamamatsu R1294, Dual MCP, QE = 4% 

White Cell, 8 passes, 40 nsec 

Hamamatsu C1370, 5 psec resolution, 
85 psec impulse response 

HP 5370A, 100 psec accuracy 

DEC LSI 11/23, Dual floppy disks 
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streak camera to eliminate any possible transient effects. The signals from 
the delay system then entered a light-tight enclosure, where the 0.355-pm 
pulse was reflected from a dichroic beam splitter, while the 0.532-pm pulse 
passed through it. The pulse in each path then passed through a bandpass 
filter for background noise rejection. They were then recombined by a 
second dichroic and were focused into the streak camera photocathode. The 
slight offset in the focussed spots due to the path splitting permitted the 
streak camera to record the two pulses simultaneously. The streak camera 
outputs were electronically read out by the video camera, and were recorded 
onto floppy disk for off-line processing. 

7 .4 Experimental Results for the CCR Array Targets 

The cube-corner reflector (CCR) array test was conducted at 
NASA-Goddard Space Flight Center during February 26 and 27, 1985. Arrays of 
four CCR's with two different orientations with respect to the optical axis 
were used. These target geometries are illustrated in Fig. 7.4. The dimen- 
sions given in the figure are accurate to within ±0.5 ram (±1.6 psec) on both 
the length and the depth of the arrays. The first array was orientated such 
that the target surface was normal to the optical axis. The second array 
was tilted so that the range spacing between adjacent CCR's was 4 mm. The 
one-way range for both target configurations was 389 meters. The RMS 
impulse response time of the streak camera receiver( 0^) was 27.7 psec. The 
laser beam divergence angles at both wavelengths were 0.1 mrad full width at 
the 10X point (0.233 mrad RMS). The nominal transmitted pulse width at each 
wavelength was 30 psec FWHM (12.7 psec RMS). 

The return pulse pairs obtained during the experiment were grouped into 
sets of 50 to 100 measurements. Typical pulse returns at 0.355- and 
0.532- pm from the normal and tilted CCR arrays are shown in Figs. 7.5 and 
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(b) Tilted Cube Corner Array 
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Fig. 7.4. Geometries of the range targets. 
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Fig. 7.5. Waveforms of Che typical pulsed signals at 0.532 and 0.355 pm 
reflected from the normal CCR array. 


7.6, respectively. With the picosecond time resolution capability of the 

streak camera receiver, the small-scale structure within a single pulse can 

be observed. This structure is due primarily to speckle induced fluctuation 

which results from the coherent interference of the pulses reflected from 

different CCR's. From Eq. (3.23) the speckle correlation length is /2 o , 

8 

which is much shorter than the range spreads of the pulse returns from the 
tilted array. Consequently, the time-resolved speckle induced fluctuation 
in these pulses are particularly severe as shown in Fig. 7.6. 

By using Eq. (7.2) and the system parameters, the mean signal at each 
wavelength for the normal array should be Gaussian shaped with an RMS width 
of 30.5 psec. The mean signal for the tilted array should consist of four 
partially overlapped pulses with their peaks separated by 27 psec. The RMS 
width of the tilted target is 42.8 psec by using Eq. (7.4). 

The mean waveforms in each data set were computed by averaging the 
measured waveforms. To compensate for the random time jitter of the 
receiver, the measured waveforms were aligned using a correlation technique. 
This was done by first selecting a reference pulse pair which had the 
highest intensity return at 0.355 ;nn. The remaining pulse pairs were 
shifted with respect to this reference pair such that the cross-correlation 
functions were maximized. An apparent mean pulse pair was computed by 
averaging the shifted waveforms. The shifted pulse pairs were then win- 
dowed, correlated, and shifted with respect to this mean and a new mean 
pulse pair was calculated. This iteration procedure was terminated when no 
further shifting was required to maximize the correlation function. 

Normally two to three iterations were required in order to determine the 


true mean waveforms. 
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TILTED CCR ARRAY 




Fig. 7.6. 


Waveforms of the typical pulsed signals at 0. 532 and 0.355 Uni 
reflected from the tilted CCR array. 
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The waveforms of the typical mean signals for the normal array are 

shown in Fig. 7.7. As predicted, these waveforms are approximately Gaussian 

- 1/2 

shaped. The widths at the e point for the 0.532 and 0.355 um are 68 and 
61 psec, respectively, which are approximately equal to twice the predicted 
RMS pulse width. Fig. 7.8 shows the typical mean signals for the tilted 
array. These pulses exhibit some sharp glints that correspond to the reflec- 
tions from CCR' s at different range. The spacings between these glints are 

- 1/2 

between 27 and 30 psec as predicted, and the width at the e point is 

102 psec for the mean return at each wavelength. 

The measured mean pulse widths on each data set are listed in the 

second and third columns of Table 7.2. For the normal array, the measured 

mean pulse widths at 0.532 um are between 26 and 36 psec, with an average 

value of 31.4 psec. The measured mean pulse widths at 0.355 um are between 

25 and 33 psec, with an average of 29.2 psec. Both these values are very 

close to the predicted value of 30.5 psec. For the tilted array, the 

measured mean pulse widths at both 0.532 um and 0.355 un are between 37 and 

53 psec, with the average widths of 44.6 psec and 44.7 psec, respectively. 

Again, both these values agree well with the predicted width of 42.7 psec. 

2 

The power spectra of the mean signal measurements ( Jg(u>) | ) were com- 
puted using a 512-point FFT. The spectra of the mean signals shown in Fig. 
7.7 for the normal array are plotted in Fig. 7.9. The theoretical predic- 
tion using Eq. (7.7) is 


2,2 2 

- u) ( 0. + rj \ 

n 


|* s (w) | 2 = <N> 2 e 


(7.17) 
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Fig. 7.7. Waveforms of the typical mean pulsed signals at 0.532 and 
0.355 ua reflected from the normal CCR array. 
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TILTED CCR ARRAY 




Fig. 7.8. Waveforms of the typical mean pulsed signals at 0.532 and 
0.355 uo reflected from the tilted CCR array. 
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TABLE 7.2. 


PULSE 

WIDTH AND BANDWIDTH MEASUREMENTS FOR THE PULSE REFLECTIONS 
FROM BOTH THE NORMAL AND TILTED CCR ARRAY 

NORMAL CCR ARRAY 


Mean pulse width (psec) 


RMS bandwidth (GHz) 

set no. 

0.355 m 

0.532 um 


0.355 Um 

0.532 um 

1 

35.6 ±1 .0 

30.5+1.0 


3.93±0.07 

3 . 59 tO .07 

2 

31.5 

29.5 


3.73 

4.20 

3 

28.5 

29.5 


4.07 

3.39 

4 

35.6 

30.5 


3.86 

3.39 

5 

31.5 

33.6 


3.39 

2.44 

6 

30.5 

25.4 


4.15 

3.12 

7 

26.4 

25.4 


3.95 

3.52 

average 

31.4 

29.2 


3.87 

3.38 



TILTED CCR 

ARRAY 




Mean pulse width (psec) 


RMS bandwidth (GHz) 

set no. 

0.355 ura 

0.532 um 


0.355 um 

0.532 um 

1 

52.9±1.0 

52.9±1.0 


2.04+0.07 

2 . 24 to . 0 7 

2 

50.8 

47.8 


2.44 

2.76 

3 

38.8 

38.5 


2.56 

2.44 

4 

40.8 

39.6 


2.04 

1.69 

5 

42.5 

43.7 


2.82 

2.32 

6 

37.1 

37.6 


2.49 

2.08 

7 

45.3 

46.8 


2.04 

1.76 

8 

49.3 

50.1 


1.82 

2.20 

average 

44.7 

44.6 


2.28 

2.19 
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Fig. 7.9. Comparison of the predicted and measured power spectra of the 
mean signals (| $ g(u> )[ for the pulse reflections from the 
normal CCR array. 
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Eq. (7.17) is also plotted using the system parameters for comparison. At 
low frequencies the measured spectra are very similar to the predicted 
spectra. At higher frequencies the noise components, which arise from 
slight misalignments in computing the mean signals, distort the spectra 
slightly. However, the noise levels are at least 30 dB lower than the 
signal spectra and can be neglected. Notice from Eq. (7.17) that the 
power spectrum is directly proportional to the magnitude of the system 
transfer function ( I H( u») | ) when the system response width (a^) is much 
larger than the transmitted pulse width (o^). Consequently, we can estimate 
| H( oj) | using the power spectrum of the mean signal when is small and when 
the target range spread is negligible. 

The power spectrum of the mean signals reflected from the tilted array 
is 

4 

l <* 

m= 1 

(7.18) 

where r is the range spacing between adjacent CCR's. Both Eq. (7.18) 
the spectrum of the mean sign* Is shown in Fig. 7.8 are plotted in Fig. 7.10 
for comparison. From this figure, we can see chat the measured results 
3j t ee in general with the theoretical prediction at low frequency. Because 
of the large speckle fluctuations in the returned pulses from this array, 
the measured signals were not perfectly aligned, which results in the high 
frequency noise components in the spectra. 

The RMS oandwidths of the mean signals can be computed using Eq. (7.9) 
and the system and target parameters. ' r he predicted values are 3.87 GHz 
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for the normal array and 2.74 GHz for the tilted array. The measured 

results at both wavelengths are shown in the last two columns of Table 7.2. 

For the normal array, the 0.532-ym mean signal bandwidth ranges from 2.4 to 

4.2 GHz with an average value of 3.38 GHz, and the 0.355-pm mean signal 

bandwidth ranges from 3.4 to 4.2 GHz with an average value of 3.87 GHz. For 

the tilted array, the 0.532-pm mean signal bandwidth ranges from 1.7 to 2.8 

GHz with an average of 2.19 GHz, the 0.355-pm mean signal bandwidth ranges 

from 1.8 to 2.8 GHz with an average of 2.28 GHz. 

2 

The mean power spectra (<|S(w)| >) in each data set were computed by 

averaging the power spectra of the measured signal waveforms. The typical 

results are plotted in Fig. 7.11 and 7.12 for the returns from the normal 

and tilted CCR array, respectively. The power spectra of the mean signals 

are also plotted for comparison. The flat region in the high frequency band 

2 

corresponds to the speckle and shot noise levels. As expected, < j S | > and 

2 2 

| are very similar at low frequency. Consequently, < j S j > can be used to 

J 

2 

estimate | $ | when the mean signals are not available. At high frequency, 

w 

however, the shot noise and speckle dominate the mean spectra so that they 

2 

are very different from j ♦<, | . 

The atmospheric delay was predicted using the spherical model for the 
atmospheric refraction and the system parameters. The detailed calculation 
is shown in Appendix C. For the round-trip horizontal path of 778 m and for 
the operating wavelengths used in the experiment, the predicted differential 
propagation time is 59.7 psec. The differential propagation time was esti- 
mated by cross-correlating the measured pulse pairs and searching for the 
time occurrence of the correlation peak. In Fig. 7.13 the estimated results 
are compar- 2 * 4 with the prediction for the two-color returns from the normal 
CCR array. The error bar corresponds to il standard deviation. During the 
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NORMAL CCR ARRAY 




Fig. 7.11. Comparison of Che measured results on the power spectra of the 
mean signals ( |* s (u>)| 2 ) and the average power spectra of the 
measured signals (<|S(cj)|^>) for the pulse reflections from 
the normal CCR array. 


POWER SPECTRUM (DB) POWER SPECTRUM (DB) 


179 


TILTED CCR ARRAY 




Fig. 7.12. Comparison of the measured results on the power spectra of the 
mean signals ( I * g (t*» ) I 2 ) a ^d the average power spectra of uhe 
measured signals :> (<j S(oj )| ^>) for the pulse reflections from 
the tilted CCR array. 
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Fig. 7.13. The measured differential propagation time between signals at 
0.532 and 0.355 pm using cross-correlation estimator for the 
pulse reflections from the normal CCR array. 
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experiment different receiver aperture sizes were used which gave rise to 
different timing accuracies. The estimated result: agree consistently with 
the predicted dispersive delay. The average bias among seven measurment sets 
is 0.48 psec, which corresponds to an average bias of 0.072 mm on the one- 
way differential pathlength. 

The estimated round-trip differential propagation tim:s for the tilted 
CCR array is shown in Fig. 7.14. Because of the broader pulse ana time- 
resolved speckle fluctuations, the accuracies of these measurements are 
worse than those of the normal target. The average bias among eight 
measurement sets is 6.1 psec, which corresponds to an average bias of 0.9 mm 
on the one-way differential pathlength. 

Fig. 7.15 compares the timing performance of the correlation receiver 

to some commonly used estimating schemes. At the time of the experiment 

there were problems with the gain sensitivity of the streak camera, and it 

was not possible to measure the signal photocounts (<1T>) accurately. 

Therefore, three different receiver aperture sizes were used in order to 

determine the dependence of the timing error on the received signal 

strength, x corresponds to the Maximum-Likelihood estimator of the form 
ML 



max arg 
T 


/ S^t) In <S 1 (t + t)> dt 


(7.19) 


- max arg 
T 


/ S 0 (t) In <S 2 (t + x) > dt 


This estimator was first derived by Bar-David [1969], and it has been shown 
to be optimal when the speckle is negligible [Tsai and Gardner, 1985]. T pK 
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Fig. 7.14. The measured differential propagation time between signals at 
0.532 and 0.355 m using cross-correlation estimator for the 
pulse reflections from the tilted CCR array. 








183 


corresponds to the peak detector of the form 


T pK - max arg [S^t)] - max arg [ S 2 ( t ) ] 


(7.20) 


Fig. 7.15(a) shows the timing errors when ranging to the normal CCR 
array. The results indicate that the timing errors are insensitive to the 
change in the received signal strength. This suggests that the system was 
operated under the speckle limited condition where the speckle SNR (K^) at 
both wavelengths are much less than <N^>. This observation is predicted 
theoretically in Chapter 5. Under these ranging conditions, the performance 

A A 

x n i s comparable to that of x , and it is much better than that of 
COR ML 

A A 

t . Since the implemenation of x requires a priori knowledge on the mean 
PK ML 

received pulse shape (refer to Eq. (7.19)), it can only be applied to the 
stationary targets. On the other hand, x is applicable to dynamic 
targets, and its performance is only slightly inferior to that using 

A 

x when the speckle effect is dominant. 

ML 

Fig. 7.15(b) compares the RMS timing errors for the returns from tne 

tilted CCR array. The reference curve (solid line) demonstrates the 

- 1/2 • ... . . . 

<N^> dependence in the shot noise limited ranging condition. It should 

be pointed out that this reference curve does not give the actual error, it 

merely indicate the magnitude between the signal strength and the speckle 

SNR. The measured results show that the signal strength in this case is 

A 

almost equal to speckle SNR. Within this region, x^ i s close to optimal 

A A 

and so its performance is much better than x^^ and x pK> It was pointed out 
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Fig. 7.15. Comparison of the measured RMS differentia! t iming ^accurac ies of 
the cro3s-correlation ( t^qr) . Maximum-Likelihood ( ^ML> » and P eak 
detection ( r PK > estimators for the pulse reflections from the nor- 
mal CCR array' and the tilted CCR array. 
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in Chapter 5, however, that the performance of can be improved by 

choosing the receiver impulse response time properly. 

7 . 5 Experimental Results for the Diffuse Target 

The diffuse target experiment was performed at NASA/Goddard-Wallops 
Flight Facility during April, 1985. The target was a flat diffuse board 
tilted 10° with respect to the incident beam. The target geometry is given 
in Fig. 7.4(c). The target dimensions are accurate to within ±5 mm on the 
width and the length, and ±1* on the incidence angle. During this experi- 
ment a more sensitive streak camera was used. The receiving system had an 
RMS impulse response time of 36 psec. The transmitted pulse width at each 
operating wavelength was 30 psec FWHM. The laser beam divergence was 2 mrad 
FW at the 10% point (0.466 mrad RMS). The one-way target range was 362 
meters . 

The speckle SNR for the signal reflected from extended diffuse targets 
is given by (Eq. (5.26)) 

K £ = A r (2 tan 0 T / X £ ) 2 (7.21) 

i 

where A^ is the area of the receiver aperture, is the laser wavelength 

and 0 T is the RMS beam divergence angle. For the system parameters used 
i 

during the experiment, the speckle SNR's at both wavelengths are on the 
order of 10^. With these large speckie SNR's, speckle noise is negligible, 
and the system performance is expected to be dominated by shot noise. 

The typical measured waveforms at both wavelengths are shown in 
Fig. 7.16. Notice that the signal intensities are relatively law compared to 
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Fig. 7.16. Waveforms of the typical pulsed signals at 0.632 and 0.355 pm 
reflected from the diffuse flat board. 
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those returns from the CCR arrays. This is a consequence of the relatively 
small lidar cross-section of the Lambertian target. Since the relative 
magnitude of the background noise increases, background radiation may intro- 
duce additional ranging error. The fluctuations within the pulse returns 

were most likely due to the receiver noise. Due to the nonlinear frequency 

* , 

mixing of the 1.064 pm and 0.532 pm pulses at the KD P (potassium dihydrogen 
phosphate) tripling crystal in producing the 0.355- pm pulses, the angular 
beamwidths of the 0.355- pm pulses are narrower than the 0.532 pm beamwidths 
[Hansch, 1977]. This difference in divergence became evident in Fig. 7.16 
because the pulse broadening is directly related to the range spread of the 
target, which in turn is related to the divergence angle. 

The typical mean return signals are shown in Fig. 7.17. The 0.355-pm 
return is approximately Gaussian shaped as predicted by Eq. (7.2). However, 
the 0.532-pm return is non-Gaussian due to the nonlinear process of the 
crystal . 

2 

It was mentioned in Chapter 7.2 that the mean spectrum (<|S(u)| >) of 

the reflected signals from diffuse targets can be used to estimate the 

2 

spectra of the mean signals ( | ♦ ( w) | ), provided both the signal strength 
( <N^> ) and the speckle SNR (KO are large. When <N^> is small, which is the 
case for the target board returns, we would expect some differences between 
the two spectra. One way to quantify the difference is to compere the band- 
widths of the two spectra. The MS bandwidth of the mean signal is given in 
Eq. (7.9). The MS width of <|S(u)| > can be defined as 

00 

/ dm <m 2 < |S(u) | 2 > 


/ dm < |s(u) | 2 > 
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Fig. 7.17. Waveforms of the typical mean oulsed signals at 0.532 and 
0.355 am reflected from the diffuse flat board. 
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After substituting Eqs. (7.3) and (7.13) into Eq.(7.22) and using the fact 
that is large, we have 

B 2 s -^r ♦ . (7.23) 

2°s 2<N>0 h 

^ 2 

By comparing Eqs. (7. 9) and (7.23), we can see that B is always larger than 

2 2 

Bg, and it approaches B<, only when <N^> is large. In this experiment <IT> 

is on the order of a few hundred. Consequently, we would expect the 

-2 2 
measured B to be larger than B . 

w 

2 2 

Typical measured results on both ! | and <|S| > are shown in Fig. 

2 

7.18 for comparison. As expected the widths of <|s| > are larger than those 
2 

of $ J at both wavelengths. The 0.355-um pulse has a larger bandwidth 

w 

because it is narrower. 

The predicted differential propagation times between the two signals 
using spherical refraction model is 56 psec. In Fig. 7.19 the measured 
results using cross-correlation algorithm are compared with the prediction. 
In this figure each point corresponds to a set of 50 to 100 timing measure- 
ments. The error bars correspond to ±1 standard deviation from the mean. 
Since the received signal strength was varied by changing the MCP setting of 
the streak camera, the RMS timing accuracies vary considerably between dif- 
ferent sets. The average bias among these six measurement sets is 4 psec, 
which corresponds to a average bias of 0.6 mm on the one-way differential 
pathlength measurements. 

Because of the large speckle SNR's, speckle degradation of the timing 
accuracy is negligible. Therefore, we would expect the system performance 
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Fig. 7.18. Comparison of Che measured results on the power spectra of the 
mean signals ( | ♦vj ( ui> | - ) and the average power spectra of the 
measured signals (< |S( w) |^>) for the pulse reflections from the 
diffuse flat board. 





to be dominated by the shot noise. Within the shot-noise limited region, 
the MS timing accuracy of the correlation receiver (Eq. (5.27)) can be 
simplified to 
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P - 

<nJT + <N^> 


t 2 j (1 - % / 20 S> 
S 3 

* 2 <N 1 ><N 2 > ojj 


<N.> « K i . (7.24) 


The first term in Eq. (7.24) is the MS error due to the first-order 
shot noise effects, and the second term is due to the second-order shot noise 
contributions. Equation (7.24) is evaluated using the system and target 
parameters and compare with the experimental results. The received signal 
strengths measured during the experiment were approximately equal at both 
wavelengths. The comparison between the theoretical and the experimental 
results is shown in Fig. 7.20. As expected, the experimental results show 
the signal dependence of the timing error. The difference between the 
prediction and the experimental results is between 4 to 15 psec. This dis- 
crepancy is primarily due to the difference in divergence angles, which 
causes the footprint intensities to be different and degrades the correlation 
between the two signals (refer to Chapter 6 ). The experimental results 
approach the predicted values as the signal strength increases. This is 
because the theoretical prediction is the lower bound of the differential 
timing error, and it was pointed out in Chapter 5 that the actual timing error 
will approach this bound when the SNR is large. 

A typical error versus correlation coefficient plot is shown in Fig. 

7.21. The vertical scale in this figure corresponds to the difference be- 
tween the estimated differential propagation time and the mean differential 
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Fig. 7.20. Comparison of the predicted and measured RMS differential 

timing accuracies of the cross-correlation estimator for the 
pulse reflections from the diffuse flat board. 
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Fig. 7.21. Error of the measured differential propagation time versus 

correlation coefficient between 0.532 and 0.355 ym signals for 
the pulse reflections from the diffuse flat board. The error 
is the difference between the mean and the actual differential 
propagation time at each measurement using the cross- 
correlation estimator. 
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propagation time in each measurement by using the correlation algorithm, and 
the horizontal scale corresponds to the correlation coefficient of the measured 
pulse pair. We can see that the error appears to be smaller when the two 

» 

signals are better correlated. Consequently, the timing accuracy may be 
improved by discarding the pulse pairs that have low correlation. 

The comparison between the timing performances of different estimating 

A 

schemes is shown in Fig. 7.22. The timing error of x is larger than that 

LUK 

A 

of t . This is expected because there is an additional second-order error 
ML 

A 

associated with T„q R due to the correlation of two noisy signals. The 
second-order error, however, can be reduced to negligible level by selecting 
the proper receiver bandwidth. The peak detector has the worst performance 
due to the presence of shot noise and speckle, which makes the peak intensity 
fluctuate randomly within the reflected pulse. 

7 .6 Conclusion 

Pulsed two-color laser ranging is an effective technique in estimating 
the atmospheric delay. However, its success relies on proper design of the 
receiver timing algorithm, short laser pulses, and high resolution 
receivers. A cross-correlation timing algorithm has been developed and its 
performance has been analyzed. In parallel with the theoretical research, 
experimental work has been underway to verify the feasiblity of this 
algorithm under actual field conditions. In this chapter, we report on the 
recent experimental timing results when ranging to the cube corner arrays 
and a diffuse flat target. Streak camera-based optical receiver was used 
during the experiment which allowed us to observe the pulse fluctuations 
caused by the speckle. For the cube corner array test, the results show 
that the differential timing accuracy improves with the increase in signal 
strength when the speckle SNR is high, and it is insensitive to the change 
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Fig. 7.22. Comparison of the measured RMS differential timing accuracies 
of the cross-correlation ( t cqr) • Maximum-Likelihood (tj^) » an d 
peak detection (t_^) estimators for the pulse reflections from 
the diffuse flat board. 




197 


in signal strength when the speckle SNR is low. These results agree in general 
with the theoretical predictions. For the diffuse target test, the results 
show that the shot noise ia the limiting factor, and the difference between 
the experimental results and the theoretical prediction is between 4 to 15 
picosecond. This difference is primarily caused by the slight difference in 
the received signal intensity at the two wavelengths. By taking this into 
account, we feel that the theoretical and experimental results have very 
good agreement. In both tests, the performance of the correlation algorithm 
was compared with that of the Maximum-Likelihood estimator. The results 
indicate that both estimators perform almost as well in the speckle limited 
region, and the ML estimator performs better in the shot noise limited 
region. As point out in Chapter 5, however, the performance of the correla- 
tion algorithm can be improved in the shot noise limited region by 
choosing the receiver bandwidth properly. Based on the experimental results 
and those of Chapter 5, and by taking into account the versatility of the 
correlation algorithm, we do feel that it is an effective and practical 
algorithm in two-color differential timing estimation. 
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8. AIRBORNE MEASUREMENTS OF SEA STATE AND ATMOSPHERIC PRESSURE 


8 . 1 Introduction 

Atmospheric pressure over the ocean can be remotely measured by per- 
forming two-color ranging to the ocean surface [Gardner, 1979]. Recording 
the ocean reflected waveforms with high temporal resolution also permits 
sea state ani wave height information to be obtained with this technique 
[Tsai and Gardner, 1982]. In order to verify the feasibility of this tech- 
nique, initial airborne two-color pressure measurements were collected in 
1983 with a photomultiplier-based system which had 500 psec time resolution 
[Abshire et al., 1984, Tsai and Gardner, 1984]. The airborne results indi- 
cated the change in atmospheric pressure between 1000 and 4000 ft flight 
altitudes, and the high correlation between the return waveforms at the two 
colors. Unfortunately, the pressure measurement accuracy was limited by the 
relatively low receiver bandwidth. In April 2 and 3 of 1985, another air- 
borne experiment wan conducted using a streak camera-based receiver with 
80-picosecond resolution in each channel [Abshire et al., 1985]. In this 
chapter we shall present the sea-state observations and pressure measure- 
ments based on the experimental results. The pressure measurement results 
will also be compared with the theoretical prediction. 

8 . 2 Theoretical Prediction 

In Chapter 7 we calculated the average power sprectrum of the reflected 
signal 


< |S(u>) | Z > 
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where H(ui) is the transfer function of the receiver, G(u») is the Fourier 
transform of the system point target response, ♦g(u) is the Fourier trans- 
form of the mean received signal, <N> is the expected signal strength and K 
is the speckle SNR. If we take the expectation of Eq. (8.1) with respect to 
the target profile £ we obtain 


< | S ( u») | 2 > * <N>|H(w)| 2 + |g(u>)| 2 + <N> 2 |G(ui) | 2 D(u>) (8.2) 


where 


D( id) 
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(8.3) 


b 2 was defined in Chapter 3 as the footprint intensity cross-section, and 
p(£p£ 2 ) i- 3 the joint probability density function of the ocean surface 
profile. We have neglected the beam curvature effects in Eq. (8.3). Since 
the wave height at any point of the ocean is due to the combined action of 
the wave motions from different directions which propagate to the obser- 
vation point, it can be argued that the ocean wave height is approximately 
Gaussian distributed. If we further assume that the sea surface is isotro- 
pic, D(w) can be express is 
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R^ is the correlation function of the ocean surface profile and is the 
profile variance. If we mcdel as a Gaussian function 

2/o T 2 

- p /2L 

R^p) • a e (8.5) 


where L is the profile correlation length and assume that the laser 
footprint is small compared to L, Eq. (8.4) can be approximated by 

D(u>) * - j (8.6) 

1 + (u>/u) ) 
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where 
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Notice that 


S 



( 8 . 8 ) 


is the RMS slope of the ocean surface profile. The average signal spectrum 


now becomes 
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<|S(u>) | 2 > 


<N> |H( »u) | 2 + <N> 2 
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I 

1 + (u>/w o ) 2 



(8.9) 


If the system parameters are known a priori, Eq. (8.9) can be used to esti- 
mate the RMS slope of the ocean by choosing u) q so that Eq. (8.9) fits the 
experimental data. 

In order to minimize the pressure error dependence on temperature and 
elevation angle uncertainties, the two-color ranging system was reconfigured 
as an altimeter at a nadir viewing angle [Gardner, 1979]. By assuming a 
spherically symmetric model for the atmospheric refraction, the atmospheric 
correction (AC) can be approximated as 

h 10“ 6 N 

AC 2 2 / 3 , 1 dh (8.10) 

J sin E 



aircraft altitude. In Eq. (8.10), we have neglected the term due to the 
difference between the geometric length of the ray and the straight-line 
path. This term is negligible for a nadir viewing altimeter operating at 
airborne altitudes. From [Marini and Murray, 1973], the group refractivity 
is approximately given by 


N > 80.343 f ( X) £ 
g T 


( 8 . 11 ) 



where P and T are the pressure and temperature, respectively, and f ( X) was 
given previously by Eq. (2.10). In an isothermal atmosphere [Wallace and 
Hobbs, 1977], the pressure will vary exponentially with altitude: 
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exp 



( 8 . 12 ) 


where 
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RT 
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MG 


(8.13) 


is the atmospheric scale height, P g is the surface pressure, M is the mole- 
cular weight of air, R is the universal gas constant, T is the surface tem- 

s 

perature in degrees of Kelvin, and G is the acceleration due to gravity. By 
substituting Eqs. (8.12) and (8.13) into Eq . (8.10), we have 
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(8.14) 


where N is the surface refractivity 
s 


P 

N =■ 80.343 f(X) rA 
s T 

s 


(8.15) 
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For E * 90*, P g = 1010 mb, and T g = 300*K, the scale height equals 
8787 m. With these parameters, the difference between the roundtrip optical 
pathlengths measured at the two wavelengths ( AR) is 


AR = AC l - AC 2 = 4.76 (m) 


1 - exp 


h 
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[f( X x ) - f ( ) ] . (8.16) 


During the April airborne experiment, the 0.532- and 0.355- van laser 
wavelengths were used. The corresponding value of AR is 
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From Eq. (2.20), the pressure sensitivity to AR is 


3P . 0.212 sin E 

Tar ' U\) - ffxj) 
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(8.18) 


Consequently, pressure change of 1 mbar caused a differential atmospheric 
delay change of 1.34 psec. 

8 . 3 System Configuration 

The block diagram of the two-color laser altimeter used during the 
experiment is shown in Fig. 8.1 with the system parameters given in Table 
8.1. The system was developed by J. B. Abshire and his colleagues at 
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The pulsed, two-color laser altimeter used during the airborne 
experiment in April, 1985. 
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TABLE 8.1. 


SYSTEM PARAMETERS FOR THE TWO-COLOR LASER ALTIMETER USED 
DURING THE AIRBORNE EXPERIMENT OF APRIL, 1985 


System Parameters 
4/85 Flight Experiment 


Laser 

: Quantel YG40 - Mode Locked ND : YAG 

nominal 30 psec pulses at 1.064, 0.532, and 0.355 urn 
0.3 mJ at 0.532 pm, 0.5 mJ at 0.355 pm 

Telescope 

: Celestjon 14 - 36 cm dia. Schmidt Cassegrain 

891 ctn effective collecting area 

Mirror 

: CVI Laser double stack dielectric 

maximum reflection at both 0.532 and 0.355 pm 

Photomultiplier 

: Hamamatsu R1294U, Dual MCP 

QE = 4% at 0.532 pm, gain = 1.3 x 106 

Optical Delay 

: White Cell, 8 passes, 40 nsec 

Streak Camera 

: Hamamatsu C1370, 5 psec resolution at sweep spreed = 2 

85 psec impulse response in system 

Readout 

: Hamamatsu C1098, Dual channel, DMA interface 

TIU 

: HP 5370A, 100 psec accuracy 

Computer 

: DEC LSI 11/23, Dual floppy disks 
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NASA-Goddard Space Flight Center. The operation of this system is essen- 
tially the same as that of the system used during the horizontal path 
experiments (refer to Chapter 7.3). The major difference is that the 
transmitter and the receiver assembly in this flight experiment was operated 
in a noncoaxial nadir viewing configuration. The transmitter and receiver 
axes were approximately parallel to ensure detectable signal returns at 
higher altitudes. 

The system was installed on board the NASA Electra aircraft. A picture 
of this aircraft is shown in Fig. 8.2. The Celestron 14 receiving telescope 
used during the experiment is shown in Fig. 8.3. To the left of the 
telescope was the transmitting mirror which directed the laser beams down to 
the ocean surface. Both the telescope and the transmitting mirror were 
directly above the nadir viewing port. Fig. 8.4 shows the receiving optics 
assembly and the streak camera. Fig. 8.5 shows the waveform recording and 
timing units. 

Two nighttime missions were flown with this system on April 2 and 3, 

1985. The missions extended from the Chesapeake Bay just west of Wallops 
Island, Virginia, over the Bay-Bridge Tunnel, and out over the Atlantic 
Ocean. In the first mission more time was devoted to system adjustments. 
Receiver triggering problems were encountered during the course of this 
flight, and as a consequence only a few usable waveforms were recorded at 
1000 ft altitude. The triggering problem was solved prior to the second 
mission, and the waveform collection process went smoothly thereafter. More 
than 1800 waveform pairs were collected at flight altitudes of 1000, 1500, and 
2000 ft. The flight route that we followed during the April 3 mission is 
shown in Fig. 8.6. 
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A picture of the Celestron 14 receiving telescope on board 
the aircraft. 
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A picture of the receiving optics assembly and the Hamamatsu 
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Fig. 8.5. A picture of the waveform recording and timing units. 
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8.4 Data Analysis - Sea State 

The measured waveform pairs in all three altitudes exhibit pulse shapes 
that are related to the various local surface profiles of the ocean waves 
within the laser footprint. 

Fig. 8.7 shows a typical return pulse pair that exhibits very sharp 
single glints. The widths of these glints are comparable to the point 
target response width. This seems to indicate that some part of the local 
ocean surface was relatively flat within the laser footprint so that return 
signals are dominated by the specular components. These sharp returns were 
obtained more frequently at lower altitudes (1000 to 1500) where the RMS 
laser footprint diameters were on the order of 30 to 43 cm, and near the bay 
where the ocean was relatively smooth. 

Fig. 8.8 shows another type of measured waveform pairs that exhibit 
double peaks. This pulse type also occurred more frequently at lower alti- 
tudes. We believe that these returns were due to the reflections from ocean 
waves that had trochoidal (or sinusoidal) local surface profiles within the 
laser footprint. The first peaks of the returns correspond to the reflec- 
tions from the wave crests and the second peaks correspond to the reflec- 
tions from the wave troughs. Since the wave heights of these local 
trochoidal waves varied slightlv, there were slight phase delays between 
reflections from various crests (or troughs). Consequently, the widths of 
the glints are slightly broader than the point target response width, and 
there are some speckle-induced pulse fluctuations. The mean wave height 
within the laser footprint can be approximated by measuring the separation 
time of the two peaks. In this particular example the separation times at 
both wavelengths are ~ 142 psec, which corresponds to the mean wave heigth 
of ~ 2 cm. This wave height value is typical for small gravity waves 
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pulse pair that exhibited single 
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Fig. 8.8. A typical ocean-reflected pulse pair that exhibited double 
dominant glints. 
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generated by light wind as a disturbing force and gravity and surre.ce ten- 
sion as restoring forces. We thus believe that this pulse type 
corresponds to the reflections from a relatively calm local ocean surface 
that consists of a few small-amplitude trochoid gravity wa-es. 

Fig. 8.9 shows yet another kind of measured pulse shape. This pulse 
type was observed more frequently at 2000 ft altitude where the laser 
footprint diameter was about 60 cm. The returns of this kind had widths 
that were in general much broader than those of the other two that w? have 
just mentioned, and they also shaped more like a Gaussian. We believe that 
these pulse returns were probably the result when ocean waves from different 
regions propagated to the laser footprint and reflected the signals by the 
scattered specular points with slopes oriented normally to the optical path. 

In order to illustrate the effects of ocean dynamics on the return 
pulse shapes, successive pulse measurements are shown at various flight line 
locations in the following two figures. 

In Fig. 8.10 we show six pulse pairs collected successively near the 
bay and at 2000 ft altitude. These measurements were collected at the rate 
of about 1 per second. We can see that there were a lot of structural 
fluctuations between successive waveforms due to reflections from different 
parts of the ocean (aircraft speed was about 18 m/sec). This figure indica- 
tes a trend that the signal waveforms changed from double peaks to a broader 
Gaussian and so on. The same trend continued throughout the entire data set 
of 100 waveform pairs. This tended to support the argument that the sea 
around the bay area consisted mainly of the trochoid waves of long wave 
periods and the superimposed car' 1 ’ary wavelets of short correlation 
lengths. Trochoid waves have sharp crests and broad troughs so that the 
reflections from the crests are weaker than those from the crests. When the 
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wave period was long, the laser footprint first illuminated the crest and 
the beginning portion of the trough which resulted in a double-glinted 
return signal. As the aircraft moved, the entire trough was illuminated 
next, which resulted in a Gaussian return. This Gaussian return was broader 
in time extent because of the range spread of the capillary wavelets and the 
possible slight tilt of the aircraft with respect to nadir. The pulse 
intensity fluctuations in successive waveform pairs were due to the combined 
effects of change in ocean dynamics, difference in laser footprint inten- 
sity, shot noise, and tilting of the aircraft. It should be noted that the 
aircraft instability (tilt) had a two-folded effect on the return signals. 
First, it attenuated the specular signal intensity; and second, it caused 
the delay (broadening) on the peak of the scattered component (refer to 
Chapter 5.3). The position shifts between successive waveforms were pri- 
marily due to the receiver time jitters. 

In Fig. 8.11 we compare the pulse shapes between measurements collected 
around the bay and those collected over the open ocean. The first three 
points were measured near point 5 and the last two were measured near point 
9 of the flight line (refer to Fig. 8.6). The altimeter was at 2000 ft 
flight altitude. The shapes of those pulse pairs collected near the bay 
followed the typical waveshapes given in Figs. 8.7 through 8.9. As the 
aircraft approached the open ocean, however, the received signal intensity 
dropped significantly and the pulse shapes changed randomly. The decrease 
in signal intensity was probably due to the increase of foams and white caps 
generated by the stronger sea wind. The white caps covered part of the 
local ocean surface, and so the number of specular points within the laser 
footprint became less. In this case, the ocean acted more like a diffuse 
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target than a mirror-like reflector. The sea wind also perturbed the corre- 
lation of the ocean surface, which resulted in more fluctuations on the 
return signals. When ranging to an open sea, special attention must be 
given to the alignment of the two laser beams in order to assure goou inten- 
sity correlation of the return signals. 

The average power spectrum for the 0.355-pm returns in each data set 
was computed by averaging the power spectra of the measured waveforms. The 
results were used to estimate the RMS slope of the ocean surface profile 
(S). To do so Eq. (8.9) was used to curve fit the measured spectrum at the 
-8 dB point by using the system parameters and adjusting the RMS slope. The 
relationship between the -8 dB bandwidth (B_g) and oj o is given by 
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A typical result is shown in Fig. 8.12. The measured spectrum in this 
figure corresponds to the 0.355-pm returns collected at 1000 ft flight alti- 
tude and between point 7 and 8 of the flight line. The low intensity com- 
ponents of the measured spectrum was due to the various noise sources. The 
estimated RMS ocean slope for this particular example is 0.192. The 
measured ocean slope can be used to estimate the average ocean wind speed by 
applying the empirical relationship derived by Cox and Munk [1954] 


S 


2 


0.003 + 0.00512W. 


( 8 . 20 ) 
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Fig . 0.12. 


Curve fit the measured average power spectrum (<|s(u)|^>) by using 
Eq. (8.9) and the system parameters. 
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W is Che average wind speed in m/sec measured at 12.5 m above Che mean sea 
level. By using Che esCimaCed value for S, W is ~ 6.6 m/sec (~ 13.2 knoCs). 

The esCimaCed RMS ocean slopes and wind speeds aC differenC parCs of 
Che ocean are ploCCed in Fig. 8.13. The resulcs were obcained aC differenC 
flighc alcicudes. The resulCs show chac Che ocean profile had a fairly 
consCanC slope of ~ 0.12 near Che bay (poinC 1 Co poinC 4), and goC larger 
Coward Che open sea. This is reasonable because Che wind over Che open sea 
is usually sCronger. From Chese values, Che esCimaCed wind speed was abouC 

4.5 knoCs around Che bay area, and ic gradually increased Co abouC 14 knoCs 
aC poinC 9. As expecCed, Che measured ocean slopes do noC appear Co be 
sensidve Co Che change in alciCude. 

AlChough we did noC evaluaCe Che sensicivicies of Che inferred ocean 
slope wich respecC Co Che various sysCem paramecers, Che inicial resulC, 
however, indicaCes Chat Che specCral Cechnique may be a possible way of 
esCimaCing Che sea sCaCe and cross-wind profile over Che ocean. 

8 . 5 DaCa Analysis - Pressure Sensing 

The cross-correlacion algorichm was used Co esCimaCe Che differencial 
propagacion Cime of Che laser pulses aC 0.355 and 0.532 pm. The measured 
resulcs were compared wich Che predicCed Cwo-color aCmospheric delay and Che 
pressure aC differenC alcicudes. 

Fig. 8.14 shows Che pressure measuremenC resulcs aC Che mean alciCude 
of 1064 fc . The lengCh of Che arrow indicaCes Che duracion of each sec, 
while Che arrow head indicaCes Che forward direcdon. The error bar 
corresponds Co ±1 normalized sCandard deviacion (o^) of the form 

a = o//n 


( 8 . 21 ) 
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Fig. 8.13. The estimated RMS sea slope (S) and the ocean wind speed (W) using 
the spectral curve fitting technique. 
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where n is the total number of measurements and a is the single-shot stan- 
dard deviation in each data set. The result shows that three of the 
measurement sets at this altitude had errors less than 3.5 psec (2.6 mbar), 
four sets had errors between 16 and 18 psec (12 and 13.5 mbar). The set 
that had the largest error of 27.5 psec (20.5 mbar) was collected during the 
time that the aircraft encountered very strong turbulence. The average bias 
with respect to the predicted differential atmospheric delay among these eight 
measurement sets was 12.4 psec (9.25 mbar). 

Fig. 8.15 shows the pressure measurement results at the mean altitude 
of 1535 ft. In these results, two of the measurement sets nad errors less 
than 2 psec (1.49 mbar), two sets had errors of 4 to 5 psec (3 to 3.73 mbar), 
and two sets had errors between 7 and 9.5 psec (5.2 and 7.1 mbar). The 
average bias among these six measurement sets was 3.58 psec (2.67 mbar). 

Fig. 8.16 shows the pressure measurement results at the mean altitude 
of 2103 ft. In these results, four sets had errors less than 2.5 psec 
(1.87 mbar), and the remaining sets had errors between 4 and 6 psec (3 and 
4.48 mbar). The average bias among these nine sets was 1.72 psec 
(1.28 mbar). 

For the pressure measurement results at all three altitudes, 10 of 23 
sets showed errors less than 4 psec, 17 of 23 sets showed errors less than 
7 psec. The large biases in some of the measurement sets were most likely 
due to the different intensity footprints of the 0.355- and 0.532-pm laser 
beams . 

The RMS differential timing accuracies (single-shot) of the cross- 
correlation estimator are compared with those of the centroid '"stimator and 
peak detector in Table 8.2. We can see that the performance of the 
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I J 

TABLE 8.2. 


COMPARISON 

OF THE SINGLE-SHOT 
FOR DIFFERENT 

DIFFERENTIAL TIMING ACCURACY 
ESTIMATORS 

II 

TT 

l! 


RMS Differential Timing Accuracy 

(psec) 

it 

Estimator 

1000 ft 

1500 ft 

2000 ft 

4 J 





n 

Cross-correlation 

25.42 

24.37 

20.48 

U 

Centroid 

62.81 

59.47 

53.51 


Peak Detection 

121.39 

112.64 

99.38 

1 1 


r 
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cross-correlation estimator *'.3 better than the other two suboptimal estima- 
tors. The peak detector is the one which is most sensitive to the shot 
noise, speckle, and ocean profile fluctuations. Consequently, its perfor- 
mance is the worst among the three estimators. 

The timing accuracies of all three estimators increase with altitude. 
This was probablv due to the effects of the non-coaxial arrangement men- 
tioned in Chapter 6. For the non-coaxial systems, the received signal 
energy is directly proportional to the overlap area between the laser 

footprint and the receiver's f ield-of-view instead of following the usual 

-2 ... 
z dependence. During this experiment the laser footprint and receiver FOV 

from the 2000 ft altitude had the largest overlap area, while those from the 

1000 ft altitude had the smallest. As a result, the 2000 ft measurements 

had the strongest signal returns while those measured at 1000 ft altitude 

had the weakest (refer to Fig. 6.11). Furthermore, it was shown in Chapter 

5 that the RMS differential timing error is inversely related to the 

received signal strength. This was why the differential timing measurement 

error decreased from 1000 ft to 2000 ft for all three estimators. 

The single shot measurement accuracy of the correlation estimator is 
between 20 and 26 psec. By averaging over 100 shots, we can improve the 
accuracy to 2 to 3 psec, which corresponds to the pressure measurement 
accuracy of less than 1.3 mbar. 

8 . 6 Conclusions 

In this chapter we analyzed the data obtained during the airborne 
two-color altimeter experiment. From these measurements, information on 
both the sea state and the atmospheric pressure was inferred. 
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At all thi •<• altitudes, the return signal waveforms exhibited either 
single glint, 1 » glints, or Gaussian shape corresponds to different 
local ocean surface profiles within the laser footprint. Pulse spreading 
occurred more frequently at the higher altitudes and over the open ocean as 
expected. The RMS slope at different parts of the ocean was estimated using 
a spectral curve fitting technique. The result indicated that the sea slope 
is smaller near the bay than over the open ocean. 

The atmospheric delay and pressure were estimated using a cross- 
correlation algorithm. In all, 10 out of 23 measurement sets showed 
pressure measurement errors less than 4 mbar, and 17 out of 23 sets showed 
errors less than 7 mbar. The large bias in some measurement sets were 
likely due to the different footprint intensities at the two wavelengths. 
However, the accuracies of these initial high resolution airborne two-color 
measurements demonstrated that this approach to pressure measurements is 
both promising and viable. 
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9. CONCLUSIONS 

In satellite laser ranging, atmospheric refraction is the most dominant 
error source because it increases the optical pathlengths from the ground to 
the orbiting satellite by over 2.4 m when the satellite is at zenith. Pulsed 
two-color laser ranging systems can be used to correct the refraction error by 
directly measuring the difference in propagation times between two optical 
pulses that are transmitted simultaneously at two frequencies. However, proper 
design of the receiver timing algorithm is necessary in order to obtain 
picosecond differential timing accuracies which are required for most 
applications. The requirement is further complicated by the various random 
noises associated with laser ranging such as shot noise and speckle. The study 
showed that speckle noise is particularly severe when ranging to retro- 
reflector equipped satellites due to the coherence interference in signal 
reflections from retro-reflectors of different ranges on the array. 

By taking advantage of the pulse correlation at the two wavelengths, a 
cross-correlation technique is used to estimate the differential propagation 

times. The performance of the cross-correlation technique is analyzed by 

% 

considering both shot noise and speckle. The accuracy of this estimator is in 
general a function of the detected signal strength, speckle signal-to-noise 
ratio, receiver bandwidth and the detected signal bandwidth. 

Pulse reflections from the diffuse targets, the satellite-based retro- 
reflector arrays, and the ocean were analyzed in order to evaluate 
quantitatively the performance of this estimator. For the flat diffuse 
targets, the timing accuracy is highly dependent on the receiver bandwidth and 
the characteristics of the time-resolved speckle. At low signal levels, it was 
shown that the timing performance of the correlation estimator is comparable to 
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that of the optimal Maximum-Likelihood (ML) estimator when the receiver 
bandwidth is chosen to match the bandwidth of the target reflected pulses. At 
high signal levels, however, time-resolved speckle places a fundamental limit 
on the performance of the correlation estimator. For the retro-reflector 
arrays, the timing performance is dominated by partially developed speckle so 
that the differential propagation time cannot be resolved to better than the 
pulse width of the detected signals. For the ocean, the timing accuracy is 
highly dependent on the local ocean surface profile and the surface correlation 
length. It was shown that picosecond timing accuracy is only feasible over 
regions where the surface correlation length is larger than the laser footprint 
radius. Based upon the analytical results, it appears that the overall 
performance of the two-color satellite laser ranging systems is dominated by 
the receiver bandwidth, transmitted laser pulse width and the pulse repetition 
rate and not the detected signal strength. 

The results for the retro-reflector arrays are of particular interest in 
this dissertation because both flat and spherical retro-reflector arrays have 
been deployed on satellites and the moon and are routinely usti as targets for 
geodetic measurements. Chapter 5 serves as the first analytical treatment on 
this target configuration by considering coherent fading. 

Experimental data analysis is also an integral part of this dissertation. 
In Chapter 7 the two-color ranging measurement made over the horizontal paths 
were analyzed in order to verify the theoretical results of Chapter 5. The 
theoretical and experimental results appeared to agree with each other in 
general. In Chapter 8 the ocean-reflected pulse measurements were analyzed in 
order to infer the sea state and to verify the pressure measurement technique. 
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The reflected pulse shapes were in general related to the local sea state 
within the laser footprint. A spectral curve fitting technique was used to 
estimate the sea slope and wind speed. The results showed an increase in both 
the sea slope and wind speed as the aircraft approached the open ocean. This 
spectral technique is a new approach to the wind speed estimation, and the 
initial inferred results appeared to be reasonable. However, the wind speed 
sensitivities to various parameters must be investigated in order tv de'ervcne 
the accuracy of this technique, and more experimental studies are needed ii 
order to verify its feasibility under a broader range of wind conditions. The 
measured results in atmospheric pressure were in close agreement wi. . the 
theoretical predictions at all three flight altitudes. Thus, it demonstrated 
that the two-color ranging technique in atmospheric pressur measurement is 


both promising and viable. 
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APPENDIX A. MEAN AND AUTOCOVARIANCE OF THE REFLECTED SIGNALS FROM A CUBE 
CORNER REFLECTOR (CCR) ARRAY 


The complex amplitude of the pulsed laser beam at the transmitting 


telescope is 


U T (£,t) * f(t) A T (r) exp (jut) 


where f(t) = laser pulse amplitude, 

a T (£) = complex amplitude cross-section of the laser beam, and 
id = laser frequency. 

If the CCR array is at a large range, the field incident on the array 
can be calculated using the Fresnel diffraction formula 


2 

U. ( p,z,t) = T f(t - — - § — ) a(p,z) exp [j(u>t - kz - ^ — p^ ) ] , 
1 — a c Zcz — Zz 


(£,z) = -77 / d 2 £ a (r) exp [-j j (-5 £•£) 


. k ( r 

J 7 ( T 


T^ is the one-way intensity transmittance of the atmosphere, k is the wave 
number, and \ is the laser wavelength. In deriving Eq. (A. 2) we have 
assumed that the RMS laser pulse width (a^) an d the area of the transmitting 
telescope (A^,) satisfy the condition 


co f » A T / 


(A. 4) 


If there are M CCR's illuminated by the incident beam, the reflected 
field at the reference plane immediately in front of the array can be 


approximated by the formula [Wang, 1984 J 


235 


U 2 (£,Z,t) 


M 

* I 

m=l 


V2 


ra 


m 


, 2?(p) 

Ui(p,z,t ^ J «(_£ ~ _P, 


)e 


il 


(A. 5) 


where $ is the random phase angle associated with reflection from the m-th 
m 

CCR. The $ ' s are assumed to be mutually independent and uniformly distri- 
in 

buted over (0,2ir). £(_p) is the range displacement of the reflecting surface 

with respect to the reference plane, and the remaining quantities in Eq. 

(A. 5) are defined in Chapter 3. The field in the plane of the receiving 
telescope is 


U R (r.,z,t) = a n^£.» z » t ^ ex P (j ut ) > 


R — 


(A. 6) 


a R (r,z,t) = exp [-j(2kz + ^ r 2 )] \ 3/ /2 A m f(t - 

m— 1 m 


V2 


(A. 7) 
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a(p ,z) e exp 
— m 


P m P.r 

m — m _ 

-jk( 2C(P ) ) 

2z z 


We assume the receiver f ield-of-view is adjusted so that all the signal 
energy which is incident on the telescope objective is focused onto the 
photodetector. Therefore, the t^tal received signal power is given by 

P(t) = Jd 2 _r w(jr) |A R (_r,z,t) | 2 (A. 8) 

where w(r) is an appropriate aperture weighting function. v(r) is equal to 
one for r inside the aperture and zero otherwise. The mean received power 
can be found by averaging the power over the random phase: 
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<P(t)> - z " l O m | a(_P (n » z ) \ |f(t - *£")) 


m=l 


(A. 9) 


2 2 
0*0 A / X 
m r m 
m 


(A. 10) 


o is the lidar cross-section of the m-th CCR on the array, A is its effec- 
ts m 

tive area, 6 is its reflectivity and A^ is the area of the receiver aper- 
m 

ture . 

The covariance of P(t) can be written as 


c p (t 1 ,t 2 ) = <p(t 1 )p(t 2 )> - <p(t 1 »<p(t 2 )> 


(A. 11) 


and can be evaluated by using Eqs. (A. 7) and (A. 8). By using the statisti- 
cal property of the uniform random variable 


(A. 12) 


where 


6 . . 

ij 



i*j 

i*j 


1 i= j=k=fl. 

< • 
0 otherwise 


Eq. (A. 11) can be expressed as 
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Cp( t ^ > t2 ) 


M M 

T aX z " 4 \ l °m °n l a( J2tn’ z) | 2 l a( 4i’ z) I* W( An*^n* z * X) 

m= 1 n=l 


m*n 
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l —to i —ci z — m z — n 


where W was given by Eq. (3.18). However, when the following condition 


Lfim ‘4.1 « 1 4< ’V 


V2 


V m,n 


(A. 14) 


is satisfied, Eq. (A. 13) can be simplified to 


Cp( c l » C 2^ 


T 3 4a r 2 ’ 4 ^ 1, °m °n l a( iW 2) I 2 I' 

m=l n=l 
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f(t. - tK p )) f (t. - iKo )) f (t, - '(Kp )) f(t, 
1 — tn i — n Z — m . 


(A. 15) 
^ o ) ) . 


Inequality (A. 14) is generally satisfied in satellite laser ranging where 

3 4 

the one-way range distance is of the order of 10 -10 kilometers. 

The signal-to-noise ratio of the speckle (K) is defined as the quotient 
between the square of the mean received energy and the energy variance 
[Gardner, Dec. 1977]. That is, 

K = [ J°*dt <P(t)>] 2 / J"dt 1 /" dt 2 C p (t 1 ,t 2 ) . (A. 16) 

— 00 — 00 — 00 

Eq. (A. 16) can be evaluated directly using Eqs. (A. 9) and (A. 13) and the 
result is given by Eq. (3.13). 
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For direct detection the mean and autocovariance of the signal at the 
receiver output can be calculated using Campbell's theorem (Papoulis, 1974] 


<S( t)> = <P(t)> * h(t) 


(A. 17) 


C S (t l ,t 2 ) = Ff ^°° dT <P(t)> h(t i “ h(t 2 * t) 


(A. 18) 
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T 2 ) 


where n is the quantum efficiency of the photodetector, and hf is the energy 
of one signal photon. 

The expected number of detected photons can be calculated usj.ng link 
equat ion : 


<N> - rr Q 0 T . 2 a r 2 ' 2 


(A. 19) 


with 


Q = / d 2 _p | a(_p, z ) | 2 /°°dt | f ( t ) 


(A. 20) 


and 
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1 a | a(p ,z)| 2 / / d 2 p | a(p,z) 
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G 


(A. 21) 
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Q is the transmitted signal energy, G is the gain of the CCR array. If the 
laser footprint has a Gaussian cross-section (Eq. (4.38)), Eqs. (A. 19) and 
(A. 21) can be expressed as 

M 

<N> = ( 2ir ) _ 1 Q T a 2 A r z" 4 tan“ 2 9 T £ o m (A. 22) 

m=l 


and 


G * T a / (2 tt z Z tan 0 ) , (A. 23) 

**. m l 

m=l 

respectively. <S(t)> and CgCt^^) can now be calculated readily by using 
Eqs. (A. 9) through (A. 21), and the final results are given by Eqs. (3.9) and 


(3.10), respectively. 
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APPENDIX B. EVALUATION OF THE BANDWIDTH OF THE OCEAN REFLECTED SIGNAL 


The bandwidth of the received signal defined previously in Chapter 5 is 
given by 


2 / du> u> 2 | 4> s ( o>) | 2 

/ dw| (J>g( w) j 2 


( B . 1 ) 


From Eqs. (3.1) and (5.19), the power spectrum of the expected signal is 
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• < exp [-j ■— [5 (£j) - C^)]]^ 


where we have assumed a Gaussian laser footprint with a radius of z tan 0 T 
at the exp(-l/2) point. The expectation is with respect to the surface 
profile. For a Gaussian surface profile with RMS surface height 0^, we have 


2aJ 


< exp [-j — [C(P,> - C(Pn)U> = exp 
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where *- 9 t * ie norma li zed autocorrelation function of the surface 

profile 5. 

We first assume L _> z tan 0 T and a quadratic autocorrelation function 
with 8 iven by 


R f(£l ‘£2^ = 1 




(B .4) 


where L is the correlation length of the surface profile. The case of 
L « z tan 0 T is treated later. 

By substituting Eqs . (B.3) and (B.4) into (B.2) and making the 
following changes of variables 


£d ‘ £1 " £2 ’ 


and 


£3 = (£1 + £2^ 2 


we have 


2 2 
- co a 
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| d>g( 


/ 4 . 4 a 

4ir z tan 0, 
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3 2 d 
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2 U.V 

g 2 . 2u> 

~ P d + 'id 


(3.5) 


-u,V 

* 1 * 8 

t + ^ 


L > z tan 8„ 


2 z tan 0 T + tan 0 T 

L 


L ^ z tan 0^. (B.6) 


By substituting Eq. (B.5) into Eq. (B.l), we obtain 


/it ( o ) erfc(y) 
g 


L ^ z tan 0 T . (B. 7) 


For the case when L << z tan 0^,, the autocorrelation function can be 


approximated by a delta function. Therefore, 
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Y 1 6 


L < z tan 
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where 


c(o g *4 °?> 


V2 
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2 z tan 9, 


L « z tan 0 T . (B.9) 


In this case, the expected bandwidth is 
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L << z tan 0^ . (B . 10) 


The combined results 


under these two conditions were given in Eqs. 


(5.46) to (5.48). 


244 


APPENDIX C. DIFFERENTIAL PROPAGATION TIME OVER A HORIZONTAL PATH 


In the spherically symmetric refraction model , the roundtrip optical 
pathlength (R q ) measured by a pulsed laser ranging system is given by 


R 

o 



dr 


i + io -6 n 

8 

sin 8 


(C . 1 ) 


r 

• . 

r* 

4 1 

It 


where N is the group refractivity , and 0 is given by Sneli's law for a 
8 

spherically stratified medium. It follows from Eq. (C.l) that the roundtrip 
differential pathlength at the two optical frequencies is 


AR =- 


r, 

R =2 / 

n J 


10“ 6 (N - N ) 
r sin 9 gj g 2 

o 


(C . 2) 


When ranging over a horizontal path and by neglecting the fluctuations 

in the horizontal refractive gradient, N will be constant with r. 

g i 


Consequently, we have 


AR 2 2 x 10" 6 (N - N ) / 1 dr a . (C.3) 

g L & 2 r Sin 9 
o 


After combining Eqs . (C.l) and (C.3), the differential optical path- 
length can be expressed as 


* • 
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AR 


10" 6 (N - N ) R 
8 1 g 2 i 

1 + 10~ 6 N 


2 10 -6 (N - N ) R [1 - 10‘ 6 N J 
g l g 2 °i g i 


(C .4) 


The corresponding roundtrip differential propagation time is 


T = 10" 6 (N - N ) T. [1 - 10~ 6 N 

gj 8 2 1 8 


(C . 5) 


where is the roundtrip propagation time of the optical pulse at wave- 
length For the wavelengths used during the experiment, 

N and N are 313 and 290 PPM, respectively [Abshire, et al., 1983] 

8 0 . 355 8 0. 532 
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